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Abstract 


’ ‘ Cooperative projects between RSI and nunerous South Dakota agencies 
have provided an Important means of Incorporating remote sensing techniques 
Into operational programs. In Fiscal Year 1980 eight projects have been 
completed; each project Is abstracted In the following, titled paragraphs. 

Detection of High Moisture Zones Near Interstate 90 

High moisture zi^es adjacent to Interstate 90 In western South 
Dakota are associated with highway warping and/or pavement breakup. 

Location of the high moisture zones Indicates possible areas where 
drainage tiles can be Installed to ease the problem. Thermal, midday 
and predawn, and color 'infrared Imagery were collected at three altitudes 
(305 , 914 , 2472 m) on t.,o dates (5/24/79 and 9/14/79) over a primary 
study area, a 20-m11e S’^gment of 1-90 near Rapid City. An 8-mtle 
segment of 1-90 near Chamberlain was also flown on 9/15/79 at 914 m; 
day and predawn thermal data and CIR photography were collected. 

Visual analysis and level slicing (thermal Imagery only) revealed that 
predawn, September 9 thermal Imagery (2 a.m. 1ST) appeared best for 
Seating high moisture zones In the Rapid City area. Color Infrared 
Imagery was of valuable assistance with the thermal imagery In the 
Rapid City site, and by itself In the Chamberlain area it was more 
useful than either midday or predawn thermal imagery for locating high 
moisture zones. Further analysis and field data are needed to verify 
the utility of the September predawn thermal data in the Rapid City 
area. 

Thermal Infrared Census of Canada Geese in South Dakota 

Current methods for censusing Canada goose concentrations on the 
malnstem reservoirs In South Dakota have sizeable errors. Significant 
errors occur because many geese are feeding away from the reservoir during 
counts and return only In the evening. The objective of this project 1$ 
to determine the feasibility of an operational Inventory utilizing aerial 
thermography collected during the night when all geese are resting on 
the water. The thermal infrared eraissivlty of central flyway geese was 
measured at ,962 ± ,017, There were no significant differences in 
emissivltles between species of Canada and Snow geese or between adults 
and Immature. The apparent temperature of Canada geese was measured at 
varying environmental temperatures. The results of this data were used 
In conjunction with test flights to determine ambient temperature 
constraints for successful aerial thermal Infrared data collection. 

Geese could be distinguished from the background of water or Ice on 
aerial thermal Infrared Imagery collected with a trimetal detector at 
altitudes less than 455 m A6L. There was a 10.98X error between estimates 
of goose nim^ers from thermal Imagery and the actual number of geese 
determined from counts made on aerial photographs. 


Dutch Elm 01 state Detection In Urban Environmon t 


i » 

Dutch elm disease Is fatal to nearly all elm trees. The disease Is 
transmitted by beetles which survive in the dead or dying trees. Removal 
of infected trees is the only effected method to control the disease. 
Economic loss, including removal of dead trees was $22 million in South 
Dakota in 1978. The city of Watertown and South Dakota Division of 
Forestry are cooperating with the Remote Sensing Institute to use remote 
sensing data for early detection and removal of infected trees before the 
disfttJse i:» spread. Aerial photography collected in August 1979 proved 
to be basically unusable because of camera malfunction. The film was 
used primarily for training in basic species identification. New 
photography was collected on June 28, 1980. The data has been reviewed 
and is suitable for objectives of the project. 

A Feasibility Study for Monitoring Effective Precipitation in South 
Dakota Using TIROS-N 

A study is in progress to evaluate the potential for using weather 
satellite (TIROS-N) data to monitor effective precipitation. TIROS-N 
imagery were compared with ground measurement of precipitation and soil 
moisture. No patterns associated with nwisture variations were found 
on the standard image products. Digitally enhanced imagery are currently 
being produced and will be analyzed. 

Open and Abandoned Dump Sites in Spink County 

There are numerous open and abandoned solid waste dump sites in 
existance which are utilized without sound environmental practices. 
Several government agencies are addressing this solid waste management 
problem. One major problem is determining how many and where they 
exist. Remote sensing techniques have successfully been used to 
identify the location of these sites and provide other resource 
information concerning pollution potential . 


Infimnct of Soil Rcfltctano on Landut Slonaturts of Cropt 


Ttn piwto-lnttrpt^ted strata ware avaluatad wtthln a study area in 
•as'ttm South Dakota as a naans of partitioning landsat data into 
ralativaly honoganaous ragtons. Eeononica, Statistics » and Cooparativas 
Sarvica proyidad tha ground sanpla schana in which 255 sagnants wara 
spraad acrass tha 6-county study araa» ground data* collected in last weak 
of July and first week of August, included field boundaries, crop type, 
and crop stage on aerial photographs (1:1 5840) » Com and stmflowar 
digital data ware extracted from an August 25, 1979 Landsat Computer 
Compatible Tape. Analysis of variance for com data indicated that 
there wara significant differences among strata and among crop stages. 
Multiple comparison tasting for com grouped strata that ware not 
statistically significant*, five to seven groups were Indicated, dependent 
upon Landsat MSS band. Sunflower data have not yet been analyzed*, 
however, expansion of the ground data provided a current areal estimate 
of sunflowers for use by several South Dakota agencies. Multivariate 
analytes are needed for a comparison between corn and sunflower results. 

Model Implementation Program - Lake Herman Watershed 

The Lake Herman watershed in southeastern South Dakota was selected 
as one of seven watersheds in the United States for involvement in the 
nationaV pilot Model In^lementatiOn Program (MIP). Remote Sensing 
Institute is cooperating with nirnierous other states and local agencies 
to provide baseline and continuous resource information for watershed 
management. An information system including land cover, soil series, 
slope, drainage and land treatment is being used to locate land in need 
of treatment, sites suitable for sediment control structures and estimate 
sediment delivery to the lake. 

Six-Mile Investigation Follow-On 


An investigation of Six-Mile Creek Watershed, located in south- 
eastern South Dakota, was conducted in FY79 to develop and apply remote 
sensing technigues to assess the geohydrology and the environment of the 
watershed to evaluate the Impact of a proposed dam and recreation site 
on Six-Mile Creek. As part of the original investigation, SCS geologists 
used thermal Imagery to select sites for drilling observation wells to 
provide information for designing seepage and cutoff controls, A follow- 
on to the original investigation was requested to document the results 
of the drilling. The drilling, which occurred in August 1979, verified 
the presence of permeable al 1 uvl urn and near-suf face groundwater in the 
vicinity of the proposed d«».^ 
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REWTE SENSING APPLICATIONS TO 
RESOURCE PROBLEMS IN SOUTH DAKOTA 


The Remote Sensing Institute has cooperated with numerous local, 
state, and federal agencies in South Dakota to effect the transfer of 
remote sensing techniques. The cooperative projects have been made 
possible through the funding of NASA, the State of South Dakota, and 
numerous cooperating agencies. The overall objectives embodied in 
the NASA Applications Program are: 1) introduce remote sensing techniques 
to local, state, or federal agencies in addressing certain resource 
problems of South Dakota; and/or 2) provide alternative methods of 
improving current practices or procedures. In FY 1980 eight projects 
have been completed. 


DETECTION OF HIGH MOISTURE ZONES 
NEAR INTERSTATE 90 


INTRODUCTION 

The deteriorating condition '' ' ■ ny highways in South Dakota is 
a very severe problem. The role ' .e state's highways in transporta- 
tion is further emphasized in 1980 due to the reduction of railv/ay 
services and concomitant increase in heavy truck traffic. While much 
of the deterioration of the highways is due to age and the economics 
of maintenance (i.e. inflation and lowered taxes), numerous areas of 
the 1-90 system in western South Dakota have had a warping or break up 
problem due to the geologic nature of this region. 

Fractures, joints, and bedding planes in the Pierre shale have 
contributed to differential weathering and ultimately to variable 
expansion-contraction of materials under and near 1-90 in western 
South Dakota. Several research and characterization studies have 
been undertaken to more fully understand the nature of these phenomena 
(Sruce and Scully, 1966; McDonald, 1970; and Hoskins and Hammerquist, 
1970). These studies have accomplished much in describing the morphology 
and genesis of Pierre shale as it relates to highway breakup and/or 
warping. 

It is generally accepted that the so-called highway bumps are the 
result of two major mechanisms: 1) differential expansion of weathered 
shale and clayey textured soils at the interface of cut-and-fill 
sections, and 2) very commonly, faults in the Pierre shale and the 
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resultant "gouge" zones of 2:1 clays. The presence of soil moisture is 
a triggering factor in bump formation. Hoskins and Hammerquist (1969) 
reported water in 4 of 9 trenches which were "... dug in a rather 
scattered manner ..."; some perched water tables were also noted. One 
source of water is the dam-like effect of bumps in highway ditches, 
i.e. drainage from highway ditches is blocked by the swelling of the 
gouge material. !n many cases, the gouge areas expand more and faster 
than the surrounding material (Harrenerquist and Hoskins, 1969). Primary 
and secondary discontinuities also provide mechanisms for water movement 
in the Pierre Formation (Bruce and Scully, 1966). 

McDonald (1970) has outlined three possible approaches to eliminate 
water accumulation in the lower portions of roadbeds located in Pierre 
shale problem areas. Very briefly the methods are: 1) an asphalt blanket 
over the in-slope, ditch bottom, and part of the backslope; 2) a drain 
or tile utilizing granular materials; and 3) increasing the gradient 
in the ditches to facilitate better drainage. 

Assuming that the above methods of combating the high moisture zones 
associated with highway bumps are successful, it follows that these 
procedures should be implemented at the earliest sign of a problem. 
Hoskins et al . (1971) reported the usefulness of low altitude color and 
color infrared (CIR) photography to locate fracture zones associa-ed 
with 1-90 near Rapid City. Other remote sensing studies have indicated 
the potential of thermal infrared imagery for locating zones of high 
moisture (Rib, 1975; Stallard and Myers, 1972). 
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The major objective of this study was to Investigate the utilization 
of CIR photography and thermal infrared imagery as tools in locating 
problem areas related to highway warping or breakup associated with the 
Pierre Formation. Other objectives included a comparison of the remote 
sensing variable (i.e. altitudes, CIR and thermal, date of collection) 
and the transfer of the data to DOT personnel . 

STUDY AREA 

A 20-mile segment of 1-90 near Rapid City, representative of 
highway breakup problems, was selected as the primary study area 
(Figure 1). Fracture zones in parts of this area have been previously 
documented (Hoskins et al., 1971). This area is in a warm, very dry 
plain with an annual rainfall of 35 to 43 cm. Soils have formed in 
Pierre shale and in high terraces. 

A secondary site west of Chamberlain was added to the project 
in mid-summer due to cooperator interest (Figure 2). This site lies 
in a warm, dry plain with a rainfall of near 43 cm. Soils have formed 
in the Pierre shale. The rolling to steeply breaking topography often 
exposes several members of the Pierre Formation. This region has a 
high landslide potential. 


PENNINGTON COUNTY STUDY AREA 
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Rain 

Racording 
Stalioni 

“ 0 ' t '!P 


s^KiOii>iAL PAGE 18 
OR POOR QUAI/ITY 



SOIL ASSOCIATION 
Terrace Soils 

18. Satontc, gently undulating 
Clay uplands 

24.Pierro-Samsil, undulating 
26. Samsil-Plerro/ rolling 
Badlands 

30. Badlands, sleep 


Figure 1. Location of 1-90 study section and the soils of surrounding 
region. Precipitation recording stations are also located. 
{See Table 2). 


CHAMBERLAIN STUDY AREA 


LYMAN COUNTY 


SOIL ASSOCIATIONS 
Uplands, Clayey 
NT" PS Opal " Promise 
Z3-NT Sansarc-Opai 

Uplands, Silty 
RJ Reliance 

Bottomlands, Loamy 
M9-HG Munjor-Haverson 







M9-HC 




Figure 2. 


Location of Chamberlain study area. Precipitation recording 
stations are the numbered positions in and surrounding the 
study site. (See Table 2). 


METHODS AND MATERIALS 


Thermal imagery and CIR photography were collected twice for the 
Rapid City site and once for the Chamberlain site. Further pertinent 
details of each flight are outlined in Table 1. All flights took place 
under cloudless skies. The 4.5 to 5.5 urn thermal detector was used in 
the initial flight due to malfunction in the 8.7 to 11.5 imi theraal 
detector. 

Neither of the sites had received any rainfall immediately prior 
to the aerial data collection. Precipitation recorded during the 
growing season is detailed by monthly totals and recording station 
(Table 2). 

Photo-interpretation of the remote sensing data was accomplished 
using visual analysis. The CIR photography in its original, film- 
positive format was systematically and stereoscopical ly interpreted. 
Landform, surficial drainage, soils, and vegetation were the major 
indicators used in the delineation of potential high moisture zones. 

Thermal imagery was level sliced to accentuate the apparent temperature 
zones of interest, i.e. cool areas during day flights and warm areas 
during night flights are zones of potentially high water tables. 

Positive prints (approximately 1:10,000 scale) of the September-collected 
thermal and CIR Imagery were delivered to DOT personnel. 

Limited field checking of the remote sensing interpretations was 
accomplished by DOT personnel. Field checks were visual; the correspondence 
of several interpreted high moisture zones on the September data (i.e. day 


Table 1. Aircraft data collected for 1980 DOT project. 
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Table 2. PrtM ipitatlon records for aredi around study sites (i'^79 data 
provided courtesy of Department of Water and Natural 
Resources, Pierre). 

Precipitation (cm) 


Station 

April 

May 

^ June. _ 

Ml. 

..Ayjjj. , 

^ Sept._^ 

Totals 

Pennington 

County 

15 

2.1 

2.9 

13.2 

8.7 

4.5 

.2 

31.6 

31 

.6 

2.9 

10.3 

13.7 

5.7 

0 

33.2 

34 

1.3 

3.3 

12.6 

10.1 

4.7 

0 

32.0 

43 

N/A 

1.3 

13.1 

12.4 

6.8 

0 

33.6 

46 

N/A 

N/A 

10.3 

11.6 

7.0 

.1 

29.0 

49 

4.4 

2.7 

9.1 

8.5 

4.9 

6,2 

35.8 

Lyman 

County 

01 

4.6 

5.0 

5.3 

13.2 

5.7 

0 

33.8 

08 

5.4 

N/A 

6.3 

13.9 

9.3 

0 

34.9 

13 

7.3 

5.8 

10.3 

N/A 

5.2 

N/A 

28.6 

17 

.13 

6.1 

6.2 

12.2 

6.9 

0 

32.2 

20 

4.2 

5,2 

14.7 

12.8 

7.7 

N/A 

44.6 

25 

5,7 

7.3 

10.7 

20.5 

6.4 

0 

50.6 


See Fig. 1 and Fig. 2 for station location. 


CIR and predawn thermal) to ground features was checked in the Rapid 
City site during June, July 1980. Drilling programs to verify shallow 
water tables have not been undertaken. Consequently, interpretations 
based on the CIR and the thermal data were not fully authenticated. 
However, DOT personnel will be utilizing the interpretations when 
drilling programs are undertaken \n the Rapid City study area. 
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RESULlii AND DISCUSSION 

In order to fulfill thu major objective of this study, three filyht 
altitudes, day and night flights, thermal and CIR data, and two flight 
periods were included, It was established, within reasonable limits, 
that the remote sensing data had merit. In accord with the major 
objective. This was accomplished by: 1) reference ts a previous project 
In the Rapid City study area (Hoskins et al., 1971); t) position of 
Interpreted areas to highway patches; and 3) communication with DOT 
personnel who were familiar with the area. The next evaluation was the 
selection of the best set of variables, i.e. altitude, etc. 

The comparison was accomplished in a two-step procedure. First, 
possible high moisture zones were located oi. the remote sensing data; 
this was relatively easy since several problem areas In the study site 
had previous documentation (Hoskins and Hammerquist, 19G9; Hoskins ot 
al., 1971), The second phase of comparison was photo-interpretation 
of various highway segments where high moisture zones were indicated. 
Qualitative comparison among the remote sensing variables took place; 
this was «rcompl1shed by answering the following two questions; 

1) was the possible high moisture zone recorded on all the Imageries, 
and 2) if It was, how easy was It to locate. The Indicators of 
possible high moisture zones were presented to DOT personnel for their 
use with photographic prints. 
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The comparison of the variables was accomplished strictly in a 
visual manner. The results for the Rapid City site are listed as 
follows, by the "best" variable: 1) altitude-914 meters j X) remote 
sensing data-predavm therr.ial; and 3) season of flight-early September. 

The selection of the most useful variable Is obviously limited by 
the number of variables. Other factors also limited the extent to 
which the foregoing resuits could be extrapolated. 

The results apply only to the variables tested (e.g. a IbOO m 
altitude may be a bettor compromise between cost and resolution than 
914 meters). The results are applicable only within the constraints 
of the antecedent climatic events (Table ?). Two thennal detectors of 
different spectral response were used, thereby, causing extrapolative 
problems In comparing the spring vs fall data. However, acc|u1s1t1on of 
thermal data in the fall has been shown to be useful for locating shallow 
aquifers (flyers and Moore, 1972). Additionally, the CIR photography, 
while not rated as the "best" variable, provided complementary 
information to the thermal data as well as a means of registering the 
thermal data to ground locations. 

Knowledge of basic soil thermal properties and characteristics Is 
useful In Interpreting the thermal data. The conductivity and storage 
of heat in the soil is differentiated to a large degree by the content 
of soil air and soil moisture. A large volume of soil air is a result 
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of sen W)fsture {this is variable with texture) and as such the 
ttermi ciHuluctivity and heat capacity are reduced. Consequently, the 
Mat capacity in a d>y soil is diminished, and the transfer of the heat 
tc Meper portions of the soil 1$ restricted, relative to a case where 
soil moisture is abundant. As a result (i.e. when con^arlng "dry" to 
soils) the naximm day soil ten^erature will occur In the dry 
soil because heat capacity and thermal conductivity are low. Whereas, 
at night when the soil is radiating to a clear sky, the wet soil will 

be warmer than the dry soil because of greater heat storage (Van Wljk 

and DeVries, 1963). In other Mrds, soil with a water table within 
capillary rise of near the surface will be cooler during the day and 

wareer at night than a soil without such a water table. It must be 

noted, however, that daytiire thereial patterns are strongly influenced 
by actively transpiring plants and the direction and steepness of 
slo(^s, amoing other considerations. 

The use of CIR photography was very useful in recording differential 
vigor of vegetatiw. In this study, we looked for non-all uvial pockets 
of actively growing vegetation to indicate a favorable moisture regime, 
or shallow water table (Fig. 3). Growing vs non-growing vegetation is 
easily distinguishable on CIR photograj^y. If seepage or a shallow 
vmter table has occurred, the vigor of vegetation drawing upon the 
arkfitional soil moisture was often easily perceived on the September 
CIR photograirtry and thenwl iraag&ry (Fig. 4), The daytime thermal 
patterns also record this cooler surface, but the thermal imagery is 
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Actively growing vegetation (A) is recorded by CIR photography 
(14 September 1979). Anomalous, nonalluvial vegetative growth 
often indicates the availability of additional soil moisture. 
Magenta or red tones are actively growing vegetation. Blue- 
green is bare soil and is of a smoother texture than rangeland 
Light-toned areas are harvested crops. Scale is approximately 
1 :5500 (18 cm/km) . 




September 14 day thermal; dark is cool, white tones >39®C 

Growing vegetation on a sideslope (A) is shown on CIR photography 
and day and predawn thermal imagery. Highway patches in this 
vicinity indicate this high moistu;*»» zone may be contributing to 
the highway problem. A small saline seep is noticeable on the 
CIR photography (B). The apparent temperature of the ground is 
represented on the thermal prints by six tones, each of which 
relates to a range of temperatures. Approximate scale is 
1:5500 (18 cm/km). 
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not limited by the necessity of having a transpiring surface. The May 
CIR photography (not shown) recorded the Influence of the spring snow 
melt on vegetative growth In the road ditches; this Is not an Indicator 
of high moisture zones. Many elen»ents were considered when the 
Interpretation was being undertaken such as landscape position, slope, 
shape of the surface (I.e. concave vs convex), and other mitigating 
features. 

Some patterns on the predawn thennal Imagery and the day CIR 
photography did not correspond (Fig. 5). While the drainage pattern 
is Indicative of shale (Fig. 5a), the thermal patterns (i.e. warm areas) 
in tnis area seemed to encircle several gravelly "caps" (Fig. 5b). 
Shallow water tables are indicated by warm apparent temperatures; 
however further grou.id Investigation, including drilling, is needed 
for full corroboration. 

The complexity of vegetated surfaces Is apparent on the CIR 
photography (Fig. 6a). However, unless the vegetation is indicative 
of the subsurface phenomenon of interest, this detail is not useful 
(ad hoc). The complexity of this surface can be reduced to series of 
discrete ranges of apparent ground temperatures (Fig. 6b), but aspect 
of slope, slope, and shape of slope (i.e. concave vs convex slope) 
present confounding thermal patterns. Further reduction of complexity 
is accomplished by using predawn thermal data (Fig. 6c). The predawn 
thermal data is less biased by the confounding factors, as listed 
above, and therefore, it empirically provides the most useful 
information. The photographic tones which correspond to 12.5 to 14.5“ C 
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14 September 



Fig. 5. 
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15 September; dark is cool, 6 to 16°C 

Potential high moisture zones are revealed on the predawn thermal 
imagery but not on the CIR photography. Warm areas (possible high 
moisture zones) excluding the highway pavement are predominantly 
north of the highway in a north-south orientation. Light tones 
are warm on the thermal imagery. Scale of CIR print is approxi- 
mately 1:5500 (18 cm/km); scale of thermal prints is approximately 
1 :4500 (22 cm/ km) . 








14 September day thennal ; dark is cool 

The Influence of aspect on road cuts (A) is shown on CIR and 
thermal data; the predawn thermal is the least sensitive to 
aspect variation in this area. A potential high moisture zone 
is indicated on the predawn thermal imagery (B). Scale of CIR 
Is approximately 1:5500. Scale of thermal prints is approximately 
1 ;5000 (20 cm/ km) . 




•JKJCiNfAL* PA('.K: fS 

u‘ i-ooK (jr.-i:rv 


Fig. 7. The complexity of the topography in the Chamberlain site 

is shown. Sharon Springs member of the Piei re Formation (A) 
outcrops near the entire study area. Salinity associated with 
water movement is shown (B). Scale is approximately 1:12000 
(8.5 cmAm). 
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apparent ground temperature on the predawn, 15 September 1979, thennal 
Imagery were Interpreted as indicative of a shallow water table in tiie 
absence of other mitigating factors, e.g. open water. 

Cl>amberJ a 1 n_ S l^e 

This region provided a landscape of much variability compared to 
the Rapid City site. While the Pierre shale remained as the major parent 
material, the exposure of numerous members of the Pierre shale and the 
steepness of the topography contributed to the surficlal complexity 
(Fig. 7). Due to these factors, the CIR photography provided the more 
useful Information In the Chamberlain site. Additionally, a temperature 
Inversion (predawn flight) may have reduced the usefulness of the thermal 
data. There were no comparisons other than CIR photography vs thermal 
Imagery because only one flight took place (Table 1). Ground verification 
was limited to visual comparison of Imagery tones to ground features. 

DOT personnel are Intimately familiar with this area and they were 
able to provide much pertinent Information about this area v/hlch was 
gathered In the construction phase. Again, no dulling programs were 
undertaken, and the Information gathered from remote sensing did not 
add substantially to what was already known. 

Major Indications of high moisture zones in this area on the CIR 
photography are as follows: 1) lush vegetation on sldeslopes (Fig. 8), 

2) salinity (Fig. 7), and 3) salt-affected vegetation. 

The location of high moisture zones In both areas, other than 
alluvial soils, is not necessailly directly related to the alleviation 
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Predawn the'nial; dark Is cool 


A comparison of CIR photography and predawn thermal imagery 
is exhibited from the Chamberlain site. Lush vegetative growth 
(A) is easily distinguishable on the CIR photography but not 
apparent on the themial imagery. Warm areas on the thermal imagery 
may or may not be associated with high moisture zones. Scale of 
CIR print is approximately 1:9000 (11 cm/km) and scale of thermal 
print is approximately 1:6000. 
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of the highway breakup. Locating high moisture zones provided only 
an Indication of the moisture supply that may provide the driving force 
behind differential expansion and contraction of the soil. Ue were 
unable to locate fracture zones; however, ground experience Indicates 
that fault zones are often difficult to observe on the ground. Further 
field Investigation by DOT personnel Is scheduled In Rapid City area, 
and location of fracture zones on the Itnagery will be emphasized. 

COtICLUSIONS 

The major objective, that Is, location of high moisture zones 
with remote sensing data, was accomplished short of drilling programs. 
Comparisons among variables Indicated that ;.rcdawn thermal data taken 
In September at 914 meters provided the best Information for the Rapid 
City site. However, the CIR photography taken In September com-lemented 
the thermal data. The relation of photographic tones on thermal prints 
with an appropriate level or range of apparent ground temperature was 
determined. 

The Chamberlain site had a comparison of CIR photography and 
thermal infrared data from only one flight. The CIP photography was 
the better remote sensing product due to the confounding thermal 
complexities In the study site, i.e. parent material differences and 
steep slopes. Photographic Indications of high soil moisture on the 
CIR Imagery were: 1) salinity, 2) spots of lush vegetation in side- 
slopes, and 3} salt-affected vegetation In the highway ditches and 
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meridians. The Indication of high soil riolsture on thermal prints was 
related to warm areas on the predawn, September Imagery. 

The rentote sensing data were useful for locating Indications of 
shallow water tables. Shallow water tables or seepage are expressed on 
landscape In several fashions as well as by a changing themial rcfime. 
These Indicators were noted In both study sites, and It remains for the 
actual determination of the shallow water table (I.e. drilling to take 
place (Appendix A). Fault zones were not located on any of the Imagery, 
but further field observations hold sone promise for this determination. 


^ 
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APPENDIX A 


Deportment of Tronsportotlon 


Division of Hlghvjuoys 

Pierre. South Dokoto S7SOI 



July 10, 1980 


Mr. Kevin Dalsted 
Remote Sensing Institute 
*>outh Dakota State University 
Brookings, South Dakota 57C07 

Re: Detection of High Moisture Zones Under Interstate Highways, 

Dear Kevin: 

As part of our 3-R Interstate Program, we will be evaluating subgrade 
problems to determine any corrective action which may be necessary to 
restore the riding quality of the surface. This may even extend to removal 
of the pavement In certain locations and reworking the subgrade. 

The first project to receive this consideration will be from Box 
Elder east 9 miles on 1-90 In Pennington County. The Thermal Imagery and 
Color IR you have supplied to us for this location will be one of the 
information sources to be examined In determining the degree and type of 
corrective work. Your assistance would be much appreciated in making any 
specific reconmendations or in the interpretation of this remotr sensing 
data In order that we might make effective use of it. 

Yours truly, 


DWA:jr 

CC: M. Buhler 


MATERIALS AND SOILS PROGRAM 




Donald W. Anderson 
Materials Engineer 
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THERMAL INFRARED CENSUS OF 
CANADA uEESE IN SOUTH DAKOTA 

By 

Robert Best and Ron Fowler 
INTRODUCTION 

Canada geese concentrate on the malnstem reservoirs of the 
Missouri River In central South Dakota, during fall migration. A 
significant portion of these Canada geese are relatively large birds 
and are hunted and shot by large numbers of hunters who consider the 
birds as trophies. Commercial or controlled goose hunting operations 
are becoming more common and Increasingly efficient at providing 
geese for hunters to shoot. The goose harvest and the potential for 
larger harvest Is Increasing. This Is cause for concern by goose 
managers In both state and federal wildlife management agencies. In 
regard to (1) extent of present harvest In relation to that which 
the goose population can withstand; (2) the proper distribution of 
harvest among areas, states and people; and (3) the population size 
In relation to potential population size or population goal based on 
available nesting habitat. 

Goose population distribution, trend and status Is determined 
by visual estimates made by trained observers In small aircraft, 
during daylight flights. The trend of population distribution Is 
toward fewer and larger concentrations due primarily to Increased 
concentration 3nd availability of food as in controlled goose 
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hunting areas. Accurate estimates of numbers of birds In these 
concentrations, and thus the entire population, are becoming more 
difficult, or Impossible, to make. Photographic spot checks of 
visual estimates have Indicated that visual estimates may be 50>75 
percent lower than the actual numbers of birds in a concentration. 
Also, the degree of error probably varies among observers thereby 
resulting In Inconsistent estimates. 

Furthermore, varyl'.g numbers of geese feed In surrounding 
agricultural areas during the ^'ay, and are not present on the water 
or Inmedlate area when counts are made. A technique for estimating 
goose numbers which relies on remotely sensed Imagery collected at 
night when all geese are resting on the reservoirs would provide 
for more accurate censusing and result In a more reliable basis for 
regulation of harvest to meet population and harvest goals. 

The development of airborne thermal Infrared sensors has 
provided a potential method for censusing warm-blooded animals whose 
behavior would allow them to be more accurately censused during 
periods of darkness. Thermal Infrared Imagery has been used, with 
varying degrees of success, for the census of white-tailed deer 
(Croon et al. 1968, McCullough et al . 1968, Graves et al. 1972); elk, 
moose and deer (Wride and Baker 1977); harp seals (Lavigne and Ronald 
1975); and polar bears (Brooks 1970). Wyatt et al (1980) recently 
statistically evaluated the use of remotely sensed thermal data for 
deer census. The effectiveness depends on the thermal characteristics 
of the animal and Its habitat and the characteristics of the scanner. 
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The objectives of this project are to: (1) develop operational 
procedures to use predawn aerial thermography to census Canada geese 
{Uranta oanadanoio) \ (2) determine which of two wavelength detectors 
(4. 5-5. 5 pm, 8.7-11.5 pm) will provide maximum apparent temperature 
contrast; (3) evaluate optimal altitude-resolution parameters for data 
collection; (4) establish an ambient temperature range at which 
procedures will be most effective; (5) determine the effect of climatic 
factors on data collection Including: temperature, cloud cover, and 
wind velocity; (6) develop an emperical relationship between the 
thermal anomalies produced by the geese and goose numbers; (7) develop 
procedures to obtain contiguous aerial thermography over water resting 
areas; and (8) measure the approximate emissivlty of goose Integument. 

PROCEDURES 

The Infrared emissivlty and radiant heat loss at varying 
environmental temperatures were measured to supplement the collection 
of aerial thermography. The exact procedures used are detailed by 
Best and Fowler (1980b) which Is attached In Appendix A. 

Aerial thermography was collected five times over goose concentrations 
on the Missouri River reservoirs In the Pierre, South Dakota area 
utilizing a Daedalus* thermal scanner (Table 1). A brief description 

*Inclusion In this report of registered trade names or trade marks does 
not constitute an endorsement by the authors or the Remote Sensing 
Institute. 
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Table 1. Thermal Infrared data collection missions. 


Date 

Time 

(CST) 

Ground Level 
Ambient 
Temperature 

Sensors 

(urn) 

Altitudes 

(AGL) 

11-16-79 

1530 hr 

1 

1 

1 

1 

O 

00 

i 

8.7-11.5 

303 m 




4.5-5. 5 

455 m 





758 m 

11-16-79 

2100 hr 

8® C 

8.7-11.5 

303 m 



(10® Inversion 


455 m 



at altitude) 


758 m 

1-30-80 

1400 hr 

-5® C 

8.7-11.5 

303 m 





455 m 

1-30-80 

2400 hr 

-13® C 

00 

• 

4J1 

303 m 

1-31-80 

1130 hr 

1 

o 

O 

8.7-11.5 

303 m 


455 m 


of the characteristics of the Daedalus thermal scanner and detectors 
was reported by Best and Fowler (1980a). Aerial photography was exposed 
concurrently with the thermal Imagery during all daylight missions. 
Daylight missions were flown to develop Interpretation procedures and 
to establish the goose density relationship. The photovoltaic detector 
(4. 5*5. 5 urn) was used during the first daylight data collection mission 
only. The trimetal detector (8.7-11.5 urn) was used on all data 
collection missions. A 10® C temperature inversion was present at 
303 m (1000 ft) AGL during the 16 November 1978 night data collection 
mission. 

Overlapping data must be collected if the area of geese is 
larger than the field of view of the scanner. Contiguous coverage 
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could eas^y be collected If sufficient moonlight was present to 
allow visual navigation. Overlapping coverage could be collected on 
dark nights by using lights at farmsteads for navigation points. A 
line can be drawn between obvious thermal anomalies in the overlap area 
on the Imagery and counts can be made on opposite sides of the line on 
overlapping Images In order to prevent duplication of counts. 

Thermal data were processed Into a photographic format which was 
scaled to base maps and aerial photographs. The thermal data was 
•level sliced* Into equal temperature increments during the processing. 
The 'level slicing* process divides the voltage signal from the scanner 
Into 6 equal voltage Increments which produced discrete gray tones on 
the photographic Imagery. Each of the gray levels represent an equal 
apparent temperature Increment. Any one or more consecutive levels 
can be further divided Into 6 more levels to Increase apparent 
temperature resolution. 

Goose counts were made on enlargements of aerial photography 
using the technique developed by Chatten (1952). Ten random counts 
were made within each region of relatively uniform goose density. 

The areas of each density region were measured with a Numonics 
electronic digital planimeter and goose totals were calculated 
directly from these figures. Chatten (1952) found that the results 
from this technique were accurate to 15% or less. 

Aerial measurements of goose concentrates on enlargements of 
aerial thermography and average goose densities determined from the 
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aerial photography were used to estimate total goose numbers from the 
thermography for purposes of comparison. An average density of geese 
weighed for areal density differences was calculated from counts made 
on aerial photographs. 


RESULTS AND DISCUSSION 

The major factors contributing to heat loss from a goose are 
conduction, convection, radiation and evaporation. The temperature 
profile near the skin Is determined solely by conduction (Birkebak et 
al. 1966). A major proportion of the sensible (non-evaporatl ve) heat 
loss fro<n the back of a Canada goose resting on water can be attributed 
to a combination of convection and radiation. The radiated component, 
which can be measured with remote sensors, should be the most significant 
because the environment acts as an Infinite heat sink. Kelly et al. 
(1954) discussed the Importance of heat loss by radiation In the energy 
balance of animals. Latent (evaporative) heat loss Is Insignificant 
at low temperatures (Salt and Zeuthen 1960). Heat loss from the bill 
was negligible because resting or sleeping geese generally place their 
bill under their wing during cold weather. 

The energy emitted from the surface of a goose depends not only 
on the temperature of the feathers, but also the emissivlty of the 
surface as expressed by the Stephan-Boltzman Law (Campbell 1977): 

R ® eoT 




W 
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where: 

R » energy emitted by non-blackbody W m 

c « emisslvity of the surface 

-8 -2 4 

0 ■ Stephan-Boltzman constant 5.67 x 10 W m K 

T « absolute temperature In "K 

It Is necessary to measure the infrared emisslvity and the apparent 
temperature of the Canada goose In order to determine If sufficient 
thermal contrast exists between the goose and the background so It can 
be dlst-ngulshed on thermal Imagery. The Infrared emisslvity and 
radiant heat loss from Canada geese were measured as part of this project 
(Best and Fowler 1900b , attached In Appendix A). The mean emisslvity 
was 0.962 i .017 and was not significantly different from other species 
of geese. Canada geese have an Insulating Integument to minimize heat 
loss and reduce the temperature differential with the environment. 
However, Best and Fowler (1980b) found that Canada geese have a 
relatively high radiant temperature relative to ambient air temperature, 
which should provide sufficient thermal contrast to distinguish geese 
from the background on thermal data. The rate of radiant heat loss 
depends on the potential or temperature difference between the goose 
and the environment. The potential Increases rapidly as ambient 
temperatures decrease. 

The apparent temperature difference between Canada geese and the 
ambient temperature Is not constant and increases with decreasing 
ambient temperatures (Figure 1). This relationship must be considered 




APPARENT TEMPERATURE DIFFERENCE (*C) 
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in flight planning to determine If the apparent temperature of the 
geese will be significantly different from the background temperature 
at a specific ambient temperature. It Is an approximately linear 
relationship and can be calculated from the following statistical 
model : 


where; 


^APD 


-0.322 ♦ 6.510 


Tapo “ apparent temperature difference between 
goose and ambient (‘’C) 

* ambient temperature (”C) 


The results of a statistical evaluation of thermal data for 
deer census by Wyatt et al. (1980) Indicate that deer can be detected 
most successfully against a snow cover background which has a 
relatively uniform radiance. The background from which the geese 
must be differentiated Is open water during fall and early winter 
and water and Ice during late winter. Both water and Ice have a very 
uniform radiance and provide an Ideal background. The emissivlty of 
water Is very close to that of a perfect blackbody so the apparent 
background temperature Is very close to the sensible temperature and 
can be measured with an Immersed thermometer. 

Background temperature used in conjunction with apparent goose 
temperature data can be used to determine optimal conditions for aerial 
thermal data collection (Figure 2). The "cross hatched" region of 
Figure 2 Indicates ambient temperature - water temperature constraints 
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AMMNT TiMmMUW CONSTtAINn FM COlUCftNO 
AttlAl TMRMOMAm FOR CmSUSMO CAFIAOA OUSI 



•25 -M -IS -K> -S 0 S 10 IS 20 25 
AMMOIT TIMViRATURI (*C) 


Figure 2. Arbient temperature/background temperature constraints 
for collecting aerial thermography for detecting geese. 

The crosshatched region represents ambient/water temperatures 
where sufficient contrast would not be present to distinguish 
geese. 
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where the thermal contrast between geese and water background would 
not be sufficient to detect geese on thermal data. These data are 
based on the assumption that a minimum apparent temperature difference 
of 5° C is necessary to distinguish geese from an open water background. 

The minimum temperature that can be resolved by *^he Daedalus thermal 
scanner used in this project is O.S** C and the spatial resolution is 
0.48 m at an altitude of 303 m. Consequently, the signal recorded by 
the scanner is an average of the apparent temperature of all objects 
within a resolution cell. The results of this study indicate that 
the apparent temperature of the goose must be at least S'* C different 
from the background in order to produce a response which is different 
than the background alone. Croon et al . (1968) reports that single 
animals smaller than a fox probably can not be detected with present 
thermal scanning equipment. The objective of this project is not to 
enumerate individual geese, but to determine the areal extent of goose 
concentrations and estimate total numbers as the product of area and 
goose density as determined from aerial photographs. 

In addition to ambient temperature, other climatic factors will 
effect the radiant heat loss for Canada geese. Moon (1974) and Moen 
and Jacobsen (1974) found that there was an increase in the radiant 
temperature on the surface of the integuments of white-tailed deer 
{Odnaoilaua virgutiamw) , mule deer {O. homionua) ^ snowshoe hare 
/vt)f:i*iiKxnua) ^ cottontail rabbit {Hylivlagust fLoriilanua) , and 
red fox {Vutpuo fulva) when wind velocities increased. They reported 
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that changes at lower wind velocities had a relatively greater effect 
than changes at higher wind velocities. Parker and Harlan (1972) and 
Moen and Jacobsen (1974) found that direct-beam solar radiation would 
result in a higher radiant temperature which may increase apparent 
temperature contrasts during daylight data collection missions. Marble 
(1967) identified cloudy days as best for detecting big game animals 
from a snow background with a thermal scanner. 

Spatial resclution of the thermal scanner is a function of the 
altitude. The size of the resolution cell increases proportionally 
with the altitude. Geese could not be resolved on thermography 
collected at altitudes above 758 m (2500 ft) AGL. Goose concentrations 
could be delineated on altitudes of 455 m (1500 ft) and 303 m (1000 ft) 
AGL. Data collected at 303 m (1000 ft) AGL had the smallest spatial 
resolution and was best suited for delineating geese, especially when 
they occurred in low densities. 

Geese could not be distinguished from the background at any 
altitude on thermography collected with the photovoltaic detector 
(4.5 to 5.5 urn). The signal to noise ratio was very low, resulting 
in data with very low temperature contrasts. This may be attributed 
to the very strong thermal absolution band due to water vapor that 
begins at about 5 um (Figure 3). There is also some absorbance by 
carbon monoxide. The effects of ozone absorbtion would be minimal on 
the trimetal detector because of the low concentrations of ozone at the 
altitude flown. 


ABSORBTION {%) 






Graves et a1 . (1972) successfully located white tailed deer 
with the 3-i> urn detector under certain conditions. Croon et al . (1968) 
and McCullough et al . (1969) recommended the use of a detector sensitive 
to the 8-14 pm spectral region because peak radiation from animals 
occurs In this range, and It coincides with an "atmospheric window". 

The tone (lightness or darkness) on the thermal Imagery Is 
relative to the apparent temperatures of scene features. The lighter 
the tone the higher the relative apparent temperature. The temperature 
range on the Imagery Is determined by the operator during data 
collection. The optimal range would encompass only the total range 
of temperatures of landscapes. The thermal data can be processed Into 
analog or digital Image formats. The analog format Is a continuous 
tone Image representing a continuous temperature range. The contrast 
of the Image ran be controlled during Image generation by a variable 
gain control. 

The thermal data can also be processed Into a digital format 
with 6 equal discreet apparent temperature Increments between the 
blackbody temperatures (Figure 4). The breaks between temperature 
Increments can more easily be delineated In the digital format. In 
this case the geese could not be distinguished on the full range 
digital data which Indicates that the apparent temperature of the 
geese Is within 1.7® C of the apparent temperature of the background 
and that both fall In the same temperature Increment. 


ORlOINAi: PACK ■ 
OC POOR QUAUTY 


a. Full range digital thermal image 
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b. Level sliced digital theriial image 


Figure 4. Illustration of digital image processing and level slicing 

of thermal data. Note the increase in temperature resolution 1 

in the level sliced data. 
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One or more of these temperature Increments can be "sliced” into 
another 6 equal levels with a subsequent Increase In temperature 
resolution (Figure 4). In Figure 4 the 9.7-11.3° C temperature level 
has been sliced with a subsequent increase in temperature resolution. 

The apparent temperature increments in Figure 4b are 0.3° C as compared 
to the 1.7° C of the full range digital data in Figure 4a. The geese 
now fall in <i temperature level different from the background and can 
be identified. The geese appear in the 10.5-10.8° C apparent temperature 
range in these data. Single levels. Isolevels, can be printed for 
easier delineation. The selection of the level containing geese is 
difficult if the location of some geese is not known. 

Very low densities of geese resting on open water could be 
distinguished on thermal imagery collected with the trimetal detector 
from 1000 ft AGL (Figure 5). The temperature of the w., i background 
was slightly above 0° C and the apparent temperature of the geese was 
approximately 3.25° C. The small groups of geese appear as light 
"dots" on the thermal imagery. Larger groups have similar tones and 
appear in larger irregular shaped areas. Most groups of geese of two 
or more in close proximity could be distinguished from the water 
background. These small groups of geese could be individually 
counted on the imagery which will provide a very accurate estimate 
in low density areas. Geese resting on the shore could not be easily 
distinguisried from the background which had a less uniform radiance. 









Aerial Photograph 


Figure 5 



Full range digital thermal Imagery 


Aerial photograph and thermal Imagery (collected at 1000 ft. 
AGL with trimetal detector) Illustrating the detection of 
low densities of geese with open water background. 
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High densities of geese were easily distinguished from water 
and ice background on aerial thermography collected at altitudes of 
1500 ft or less (Figure 6). The thermal anomaly representing the 
geese has the lightest tones and warmest apparent temperature on the 
imagery. Differences in image texture within the area of oiese are 
a result of differences in goose density. Densities of geese less 
than 2000 per hectare could not be consistently interpreted on the 
imagery with snow and ice background (Figure 7). In six goose 
concentration areas on two different days only one had a significantly 
large area of low density geese. There was an 8.9X error between the 
areas of geese over 2000/hectare on the thermography and the area of 
geese over 2000/hectare measured on the aerial photography. 

The average density of geese in six goose concentrations on two 
different days, as determined from aerial photography, was 4301 geese 
per hectare. This figure was used to calculate estimates of total 
geese from area measurement on thermography. There was a 10.98% 
error difference between the estimates from thermography and the 
actual number of geese determined from counts of the aerial photography. 
Improvements in this estimate may be possible if apparent goose 
density differences were delineated on the thermography and the empirical 
goose density calculations were refined. 





Aerial photograph 


Of PQOH 


■ jfl 


High gain analog thermal imagery 


Fioure 6. Aerial photography and thermal imagery (collected at 1000 ft. 

AGL with trimetal detector) illustrating the detection of high 
densities of geese with a predominantly ice and snow background 



k* 3.086 


3,975 


3,599 


Diagram illustrating the comparison of delineation of goose 
concentration on aerial photography and thermal imagery. 
Solid line is delineation from aerial photography and dashed 
line delineation from thermal imagery. 
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SUI4MARY AMD CONCLUSIOf^S 

Canada geese can be distinguished from a water or Ice background on 
aerial thermal Infrared Imagery. The techniques developed In this project 
could be used In an operational predawn census of Canada geese during 
fall migration when geese are concentrated on refuges. Low concentrations 
of geese can be counted on the imagery. Estimates of geese In high 
concentrations can be calculated from measurements of the area of geese 
made on thermal images and empirical goose densities derived from aerial 
photographs. Level slicing of the thermal data prior to generating an 
Image Increases the temperature resolution of the Image which will 
Improve the apparent temperature contrast on digital Images. High 
gain analog or full range digital Images are the easiest to produce and 
should be used except when contrast Is very low. Computerized digital 
analysis of the data would further Increase relative contrast In the 
thermal data and Improve spatial resolution lost In the generation cf 
Imagery. 

The spatial resolution of the thermal scanner that was used In 
this project limited the altitude for data collection to 1500 ft. A6L 
or less. The best results were obtained when data were collected at 
1000 ft. AGL. The photo voltaic detector did not provide sufficient 
apparent temperature contrasts to distinguish geese from the background 
under the environmental conditions experienced during the test flight. 

Data collected with the trimetal detector produced Imagery be«:t 
suited to the objectives of this project. 
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An Infrared emIssivUy .962 was measured for Canada geese. There 
were no statistically significant differences In cmissivlty between 
species of other central flyway geese. The relationship of ambient 
temperature and the apparent temperatures of Canada geese were analyzed. 
These data were used In conjunction with cmissivlty data and the results 
of test flights to determine the ambient temperature constraints for 
successful data collection. 

These techniques could replace current census methods and provide 
more reliable census data. Accurate census data Is necessary to formulate 
harvest regulations which provide both maximum recreational opportunity 
and more equitable distribution of harvest. 
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appendix a 


1 


Robert G. Best 
Remote Sensing Institute 
South Dakota State University 
Brookings, South Dakota 57007 

INFRARED EMISSIVITY AND RADIANT HEAT LOSS FROM CANADA GEESE 

R.6, BEST, Remote Sensing Institute, South Dakota State 
University, Brookings, SD 57007 
R. FOWLER, South Dakota Department of Game, Fish and Parks, 
Pierre, SD 57501 

Ahivtraat.: The thermal infrared emissivities of 2 subspecies of 

Canada geese and 2 color phases of snow geese were measured. 

The mean emissivity of the goose integument was .962 * .017 and 
there were no statistically significant differences between the 
subspecies of Canada and the snow and blue color phases of snow 
geese. The apparent temperature of 6 hand-reared Canada geese 
was measured at ambient tenperatures varying from -25 to 20°C. 
There were no statistically significant apparent temperature 
differences between geese. All measurements were included In an 
analysis of the relationship of apparent goose temperature with 
varying environmental temperature. There was an approximately 
linear relationship (r » .87) between apparent goose 
temperature and varying ambient temperature with an increase in 
the apparent temperature difference between goose and ambient 
temperature as the ambient temperature decreased. Radiant heat 
loss was calculated from apparent temperature measurements. 
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Thcniiol Infrared emissivity and radiant heat lois 
nmasurements from Canada geese {Hmntu ihinah >u'ii<) are an 
Integral part of energy balance calculations and tan be u'ied to 
determine the feasibility for aerial thermal Infrared census. 
Radiant heat loss ineasurciiierts from Canada geese can be used to 
determine If sufficient apparent temperature contrasts exist In 
order to discern them on thermal Imagery (thermography) . Heat 
loss of 2 subspecies of Canada geese at varying environmental 
temperatures and the effect of heat loss on their distribution 
has been studied by Birkebak et al (1966) and LeFebvre and 
Raveling (1967), respectively. 

The major factors contributing to heat loss from a Canada 
goose are conduction, convection, radiation and evaporation. 

The temperature profile near the skin Is determined solely by 
conduction (Birkebak et al. 1966). The conduction Is along the 
fibers of the feathers or through air and water vapor trapped 
within the feathers. A major proportion of the heat loss from 
the surface of the Integument can be attributed to a combination 
of convection and radiation. Latent (evaporative) heat loss Is 
Insignificant at low temperatures (Salt and Zeuthen 1960). Heat 
loss from the bill will be negligible for resting geese which 
generally place their bill under their wing during cold weather. 
Only the radiated component of heat loss can be detected with 
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thermal infrared sensors. However, the radiated component can 
be very significant with the environment acting as an infinite 
heat sink. 

The thermal infrared energy emitted from the surface of a 
goose depends not only on the temperature of the feathers, but 
also the emissivity of the surface as expressed by the Stephan- 
Boltzman Law (Campbell 1977): 

R = eoT'* [1 3 

where: 

-2 

R « energy emitted by non-blackbody W m 

c = emissivity of the surface 

-8 -2 4 

a = Stephan-Boltzman constant 5.67 x 10 W m K 

T « absolute temperature in ”K 

Most surfaces radiate slightly less than a perfect radiator 
(blackbody). Hammel (1956) measured the emissivity of the 
integument of 10 species of artic birds and mammals and found 
that within experimental error, the emissivities were not 
measurably different than 1.0. Hammel concluded that there was 
no advantage In heat conservation for animals with white 
Integument. It was shown that animals which were white in 
visible wavelengths were not necessarily "white" in the thermal 
infrared wavelengths. Svihia (1956) found similar results in an 
experiment with dyed rats. 
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The objectives of this project were (1) to measure the 
emissivity of Canada geese, (2) to measure the apparent 
temperature (temperature measured remotely with a thermal 
infrared sensor) of Canada geese at varying ambient temperatures, 
(3) to determine the radiant heat loss from Canada geese at 
varying ambient temperatures, (4) to evaluate if sufficient 
apparent temperature differences are present to image Canada 
geese with an airborne thermal infrared sensor. 

The authors would like to extend their sincere appreciation 
to Mr. Gerald Jasmer for the many hours spent helping on the 
project. Wt re particularly grateful to Mr. Jean Mehegan for 
allowing us to make apparent temperature measurements of his 
hand-reared Canada geese. We are indebted to Mr. John Link, Mr. 
Paul Bultsma, Mr. Gerald Jasmer, Mr. Robert Houck and the Lake 
Sharpe Goose Ranch for donating the goose skins used for 
emissivity measurements. Funding for this project is provided 
by the South Dakota Department of Game, Fish and Parks and the 
NASA Technology Transfer Division under grant #NGL-42-003-007. 
METHODS 

Emissivity values of unpreserved goose integuments from 2 
subspecies of Canada geese aamdensin intex>iov and B.o, 

huUihinaonii) aid adult and immature snow lAfuun' 
aaevulemaens) were calculated from measurements tnade with a 
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Barnes Precision Radiation Thermometer* (PRT). The filtered 
sensitivity of the radiometer was 8-14 iim. Skins were removed 
from freshly killed geese and frozen until measurements were 
made. Skins were removed from the freezer to allow their 
sensible temperature to equilibrate with environmental 
temperatures. All measurements were made in a shaded area under 
clear skies. 

Two methods were used to estimate the emissivity. A rough 
estimate was made by ratioing the 4th power of the apparent 
temperature (°K) with the 4th power of the actual temperature 
(®K). A second more precise technique developed by Fuchs and 
Tanner (1966) for infrared thermometry of vegetation is based 
on a form of the Stephan-Boltzman Law which considers background 
radiance effects. 

R|^ = caT^ + (l-t;)B* [2] 

where: 

R|^ - radiation received by the sensor 
(l-e)B* - incoming thermal infrared reflected by the surface 

B* = background radiance 


♦Inclusion in this report of registered trade names or 
trademarks does not constitute an endorsement by the authors or 
the Remote Sensing Institute. 
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In this method the apparent temperature of an aluminum 
plate of known eniissivity is measured with the PRT. The actual 
temperature is measured with a digital thermometer and 
thermocouples entf)edded in the plate. The background radiation 
is then calculated by solving equation [2] for background 
radiance (B*). 


B* 


-j:oT 

*1-T 
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where: 

Rj^ = oT^ (T is apparent temperature measured by PRT in °K) 
Measurements of the apparent temperature of the goose skins 
were made with the PRT and the actual temperature was measured 
with the PRT by placing the goose skin under a specimen chamber 
lined with aluminum foil, which simulates a blackbody cavity. 

The emissivity of the goose skins were calculated by solving 
equation [2] for e. 


c = 


Rb - 

oT"^- B* 


[ 4 ] 


where: 


T = actual temperature °K measured with PRT under 
blackbody cavity 

Five measurements of apparent and actual temperature were 
made on each of two different skins for each type of goose. An 


I- Jv 
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average of background radiance, calculated from measurements 
made at regular Intervals during the measurement of goose skins, 
was used in the calculations. Statistical means and standard 
deviations of emissivities were calculated for each individual 
skin and goose type. Analysis of variance was used to determine 
if statistically significant emissivity differences occurred 
between geese. Statistical analyses were performed using the 
Statistical Analysis System (SAS) procedures (Barr et al . 1976). 

The apparent temperature of 6 hand-reared pinioned Canada 
geese was measured with a radiometer (PRT). Five replicate 
measurements of each goose were made at environmental 
temperatures ranging from -25 to 10°C. All measurements were 
made under cool clear skies. The background radiation was 
estimated from PRT measurements of an aluminum plate with a 
known emissivity in which thermocouples have been embedded to 
determine actual temperature. Mean apparent temperatures of 
each of the 6 geese were plotted against ambient 
temperatures. Linear regression analysis was used to 
statistically analyze the relationship. A similar technique was 
used by ftoen (1968) to measure the surface temperature and 
radiant heat loss from white-tailed deer. 

RESULTS AND DISCUSSION 

The mean emissivities of 6 types of geese ranged from .957 
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to .966 {Table 1). Statistical analysis of variance indicated 
that there was no statistically significant ernissivity difference 
between geese (P>.10). Based on these results the mean 
einissivity for all geese is .962 + .017 (95% confidence interval). 
Several authors (Hamilton, 1939 and Hesse et al . 1951) have 
hypothesized that arctic nammals and birds may have a reduced 
thermal emittance as a result of their white coloration. The 
results of this study contradict this hypothesis. Adult snow 
geese which are pure white have the same infrared eniissivity as 
the very dark colored Canada geese and both are very close to 
that of a blackbody radiator. There would be an advantage, had 
the emissivity been very low, because radiant heat loss is 
proportional to emissivity. Emissivity values calculated 
without consideration of background radiation ranged from .98 
to .99 which would correspond very closely with emissivity 
measurements of arctic mammals and willow ptormigan {lagopus 
lagopm) made by Hammel (1956). The effects of background 
radiation are not considered in the emissivity measurement 
technique developed by Hardy (1934) which was used by Hamtiiel . 

Radiant energy loss might still be small for geese, even 
with an emissivity near 1.0 if the surface temperature of the 
integument and the environment are nearly the same. However, 
the apparent temperature of a clear sky will usually be 
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Table 1. Infrared emissivitles of central flyway geese. 


Species 

Visible 

Color 

Mean 

Emissivlty 

Canada goose 

dark grey 

.962 

{BxHxnta eanadenoic intai'iox') 



Canada goose 

dark grey 

.966 

{B, hntidiinoonii) 



Adult snow goose 

white 

.963 

{Aruuxv oaavuleaaeno ) 



Immature snow goose 

light grey 

.959 

Adult blue goose 

grey 

.964 

{Anoev aaovuleooam) 



Immature blue goose 

grey 

.957 


Grand mean 

.962 


» 


» 
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significantly colder than the air temperature. Hardy and Stoll 
(1954) found that the radiant temperature of the clear sky may 
be 30 to 40“C colder than air temperature. Clear sky will act 
as a heat sink for the exposed portion of the goose. 

The relationship of the apparent temperature of Canada geese 
to varying environmental temperature is illustrated in Figure 1. 
The apparent goose temperatures were measured remotely with a 
radiometer and are an average of 5 readings from the back of 
each goose. Data for all 6 geese are Included together In 
Figures 1 and 2 because an analysis of variance Indicated that 
there were no significant differences (P>0.05) between geese for 
the apparent goose temperature/air temperature relationship. 
Surface temperatures of the geese can be determined from these 
data by using the Stephan-Boltzman law with the measured 
emlssivity (.962). The surface temperature will be very close 
to the apparent temperature as measured with PRT because the 
emissivity Is very close to 1.0 and backgrouad radiance is 
negligible. There appears to be a linear relationship between 
the apparent temperatures of the geese and the ambient 
temperatures based on these data. An estimate of the apparent 
temperature of the geese can be calculated from the following 
regression equation; 



.692 + 6.103 
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where*. 

T » apparent goose temperature "C 

Tgmb ■ ambient temperature 

2 

This equation accounts for 87% of the variance (r = .873) in 
the data. 

The difference between the apparent temperature of the goose 
and the ambient temperature is, among other variables, an 
indication of the thermal conductivity of the integument when 
the emissivity is near 1.0. These data show that as ambient 
temperature decreases the apparent temperature difference of the 
goose increases. Convective heat losses should increase as the 
temperature difference increases. 

These data can be used to determine if sufficient apparent 
temperature contrasts are present to discern geese from 
background and are essential for planning the collection of 
aerial thermography designed to image geese. 

The radiant heat loss of Canada geese 't varying ambient 
temperatures is illustrated in Figure 2. Heat loss by radiation 
can be calculated from the product of the 4th power of the 
apparent temperature (®K) as measured with PRT and 

•*8 “2 

Stephan-Bol tzman constant in the form: o = 4.93 x 10“ kcal m“ 
h“^ K“^. Radiant heat loss decreases as the apparent temperature 
of the gnose decreases with ambient temperature. It follows that 
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as ambient temperature decreases, heat loss by convection should 
Increase while radiant heat loss will decrease. 

Numerous factors can have an effect on the "apparent" 
temperature and radiant heat loss from an animal The net 
radiation loss can be affected by not only the downward flux 
under different sky conditions (Moen 1968) but can also be 
affected by absorbed solar radiation (Parker and Driscoll 1972) 
and changes In wind velocity (Blaxter and Wainman 1964, Moen 
1974, Moen and Jacobsen 1974). Knowledge of the radiation 
profile or "effective surface" for the goose Is essential for 
calculating total radiant heat loss. The radiation profile Is 
used to adjust for geometrical consideration and Is based on the 
proportion of the body exposed to different environmental energy 
flux. The radiation profile for Canada geese Is not known. 

When the radiation profile Is determined the data can be used 
In thermal energy exchange calculations. These data suggest 
that sufficient apparent temperature contrasts may be present 
In order to utilize aerial thermography for censusing Canada 
geese when environmental conditions are suitable. 
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APPLICATION OF REWTE SENSING TECHNIQUES FOR DETECTING 
DUTCH ELM DISEASE IN URBAN ENVIRONMENTS 
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INTRODUCTION 

Dutch elm disease is a disease fatal to im elm species but 
especially the American elm {Uimui'> m;rUhv,t). The disease is caused 
by a fungus {Cdwtinnjahio ulm"-) which invades and grows in the water 
conducting vessels of elm. The disease is generally transmitted by 
the European elm bark beetle (.ovoTyf-m; multi,u.inatut)) and the native 
elm bark beetle {iJylimjophnio vufip> •.•:). The beetles overwinter as 
larvae and adults under the bark of dead or dying elms. In the spring 
the emerging adult beetles, contaminated with fungus spores, fly to 
dead or dying trees to reproduce and then on to healthy trees to feed. 
It is at this point that the infection is spread. 

Sanitation by removal and destruction of dead and dying elm is 
the only way to slow disease spread within an elm population. Prompt 
identification and destruction can reduce losses to levels where 
severe economic impact is minimized. The Department of Game, Fish 
and Parks estimated that 28,000 elms were lost to Dutch elm disease 
in 1978. The economic loss, including removal of dead trees exceeded 
$22 million. 

The severity of the tree loss problem prompted an investigation 
to test the feasibility of using ren»te sensing techniques for early 
detection o^ Dutch elm disease. French and fieyer (1977) used low 
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altitude color infrared aerial photography to detect diseased trees. 
Waltz (1969) and La Perriere and Howard (1971) also recognized the 
utility of color infrared photography for detection of diseased trees. 
The investigators determined the disease is easiest to detect and has 
the highest incidence in early July. It is also recommended that a 
survey in late August would also be useful in detecting trees infected 
by second generation adult beetles. 

The concensus of the investigators was that aerial photography 
was a fast, low cost survey method. Accuracy of detection ranged from 
50-70 percent. Accuracy was influenced by difficulty of species 
identification, difficulty in discriminating the target disease from 
other kinds of moisture stress-inducing situations, difficulty in 
mapping locations of infected trees and coordinating the photo with 
the ground situation. It is expected accuracy can be improved in the 
proposed project because many of the previous investigations were 
performed in a natural forest rather than an urban forest. In a normal 
urban environment it should be easier to identify species because the 
crown density will be less. 

The city of Watertown with the assistance of the South Dakota 
Department of Game, Fish and Parks, Division of Forestry, has taken 
significant steps that should facilitate Dutch elm disease detection. 

In the past the Watertown city forester, with the assistance of the 
South Dakota Division of Forestry, has mapped and inventoried by 
species all trees on most public land and street boulevards. The 
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inventory was performed to simplify the annual ground surveys necessary 
to identify diseased trees on public land. 

The problem now facing the city of Watertown is how to identify 
early diseased elms on private property as well as to monitor for 
disease in the already identified elms on public property. The city 
has passed an ordinance permitting the city forester to condemn and 
remove diseased trees on private property at the owner's expense. 

The process of identifying the diseased trees on private and 
public land by conventional ground survey is time consuming. The city 
officials and South Dakota Division of Forestry are hopeful that remote 
sensing techniques can significantly reduce the amount of time and expense 
needed to perform the inventory of private trees. 

PROCEDURES 

It was originally proposed to collect 1:4,500 scale color 
infrared aerial photography and thermal infrared data during early 
July and again in late August, 1979. However, in 1979 eastern South 
Dakota experienced an unusually wet summer. Consequently, cloud cover 
restricted any data collection untr. late July. An aborted attempt 
was made to collect data in late July but the airplane developed 
engine problems shortly after takeoff and was subsequently grounded 
for four weeks for repairs. Color infrared aerial photography was 
finally collected in early September. A camera malfunction severely 
overexposed approximately half of each frame of the film resulting 
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In the film being of minimal use for an operational disease detection 
program. No thermal infrared was collected because of a scanner 
malfunction. 


RESULTS 

Despite the Initial problems in obtaining data both the city of 
Watertown and the South Dakota Division of Forestry were still 
interested in pursuing the project. It was mutually decided to utilize 
the data which had been collected to develop and refine interpretation 
procedures to be used in an operational program. 

A meeting was held with the Watertown city forester and a 
representative of the Division of Forestry to determine what infonnation 
and assistance was needed to utilize the existing data effectively. 

It was determined that the major use of the data would be to assist 
in training an interpreter to identify tree species. The 1979 film 
was reviewed to determine which portions of the city had suitable 
photographic coverage. The areas were delineated on a city base map. 

The map was given to the city of Watertown forester and he provided the 
corresponding inventory data.- 

The training was by compa>'ison of the photography to the existing 
public land tree inventory and other interpretation keys developed 
for tree speci^e identification. The training enabled the interpreter 
to readily identify elm trees. However, difficulties still exist 
in determining individual trees in an overcrowded stand. Overlapping 
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tree canopies is the major reason for the problems. The anreed upon 
method of dealing with the problem is to merely identify that there is 
a potential disease problem associated with a group of trees. Ground 
inspection should identify the individual problem trees, 

RECENT DEVELOPMENTS 

On June 28, 1980, 1:4,500 scale color infrared stereo aerial 
photography and thermal infrared imagery was collected for the city 
of Watertown. The photography has been reviewed and appears to be 
excellent data for identifying elm trees infected with DED. The 
thermal imagery has not been reviewed at this time. 

On June 23 a preliminary ground inspection was conducted to 
obtain ground truth information. At that time nearly 50 trees were 
identified as having DED. 

Interpretation of the 1980 photography has begun. Pink and light 
gray areas in ‘elm trees are suspected to be early indication of 
disease presence. All leafless trees are being identified as dead 
elm trees. The dead trees could be of a different species but it is 
impossible to discern this from the photography. Final determination 
will be made during the ground survey. 

The major interpretation problem is the ability to distinguish 
the light pink areas of potentially infected elms from the pink tone 
of an ash tree. The similarity is confusing but can be dealt with by 
making an accurate interpretation of species. Figure 1 is a contact 
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Fig. 1. Contact print (1:4500) of the 1980 photography with a sample 
of suspected elm trees annotated. 
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print of the 1980 photography. A sample of the suspected trees have 
been annotated. Many more will be identified in a detailed interpa*ta- 
tion using a 2x enlargement and stereoscopic viewing of the film 
transparencies . 


ANTICIPATED WORK 

Work has already begun on interpreting the photography for OED. 

The plan is to annotate DED trees on contact prints of the photography. 
The photography will also have streets and avenues annotated to enable 
ground crews to find trees suspected of having the disease. 

Currently the city forester is conducting the annual DED survey 
of public and for the first time p>'ivate land. Records of infected 
and removed trees are being kept so a comparison of the ground inspection 
and remote sensing techniques can be made. 

A followup ground inspection and aerial data collection will be 
conducted in late August. This survey will identify newly infected 
trees. More importantly, from the remote sensing standpoint, the late 
summer survey will assist in the evaluation of the potential of using 
color infrared photography for determination of previsual stress 
conditions due to DED. It is likely that some elm trees 
identified on the photography as potentially infected trees may not 
show symptoms until after the first ground inspection has been 
completed. The followup inspection will verify the usefulness of 
the data for previsual stress detection. 
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The thermal infrared imacjery obtained will be analysed to 
determine if there is any potential use of the data for discriminatirnj 
infected trees under stress. 
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A FEASIBILITY STUDY FOR 
MONITORING EFFECTIVE PRECIPITATION 
IN SOUTH DAKOTA USING TIROS-N 


INTRODUCTION 

Improved estimates of soil moisture and rainfall distribution 
are Important for agricultural and hydrological applications (yield 
forecasting, pest management, runoff modeling, flood forecasting, etc.). 
In drafting South Dakota's response to the National Climate Act, a 
Workshop of Climatic Data Users set establishment of a soil moisture 
information network as a top-priority for the state. The Workshop 
recognized that soil moisture budgeting will be an Integral part of a 
soil moisture network. 

Suitable soil moisture budget models have been developed. However, 
a weak link in the models is uncertainty in estimates of effective 
precipitation (infiltration), a major input of the models. Spatial 
evaluations of precipitation and soil moisture using ground measurements 
are difficult because of the requirement of dense networks. The 
possibility exists, however, of using satellite visible, near IR, and 
thermal IR data to monitor precipitation and soil moisture patterns 
associated with changes in evapotranspiration and thermal inertia 
as inputs into soil moisture models. 

A study was proposed to evaluate the feasibility of using weather 
satellite (TIROS-N) data to monitor effective precipitation for use 
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In a soil moisture network In South Dakota. TIROS-N collects data 
(visible, near and thermal IR) twice daily (mid-afternoon and early 
morning) at a spatial resolution of 1.1 km. 

MATERIALS AND METHODS 

Dally rainfall estimates for July and August, 1979, were provided 
by the South Dakota Department of Water and Natural Resources (DUNR). 

DWNR operates a network of approximately 1500 rain gauges. A network 
of 81 soil sampling sites was established with sites selected 
to represent the major soil associations In the state. Land use categories 
were restricted to pasture and small grains. Soil samples for determining 
water content were collected at ten day intervals for each site during 
July and August, 1979. Percent cover was also estimated at each site. 

TIROS-N Imagery were r(»co1ved and wore compared with rainfall 
and soil moisture estimates to locate apparent reflectance/omittancf? 
anomalies associated with moisture distribution patterns. 

RESULTS 

Standard TIROS-N imagery is produced to enhance cloud patterns 
at the expense of land surface features. This is particularly evident 
in thermal Infrared. As a result, no discernible patterns associated 
with moisture variations were found on the standard products. Photo- 
graphically enhanced Imagery showed little improvement over the standard 
products. 
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Selected digital data were ordered, and are currently being used 
to produce digitally enhanced Imagery. Because of delays in obtaining 
the data, and extensive reformatting required to make the data 
compatible with the RSI analysis system, the enhancement process is 
not completed. Thus, there are no results to report at this time. 

A detailed final report of results will be submitted upon completion 
of the the analysis of the digitally enhanced Imagery. 


SURVEY AND Af^ALYSIS OF POTENTIAL I'OLI.UTION 
FROM OPEN AND ABANOONFl) SOLID WASTE DUMP 
SITES USING P.EMOTE SENSING TECHNIQUES 

INTRODUCTION 

There are numerous unauthori iied open and abandoned 'jolid waste 
dumps 5tin In existence In South Dakota. The majority of these dumps 
have been and are still utilized without concern for sound environmental 
practices. Potential health and environment hazards that are created 
by open and abandoned dump sites may Include pollution of ground and 
surface waters from leachates and runoff; the breeding of germs, 
insects, and disease carrying rodents; potential fire hazard, generation 
of methane gas; and the creation of an unsightly nuisance. 

Several local, state and federal government agencies are 
cooperatively involved in solid waste managoiiiGnt in South Dakota and 
in particular Spink County. The Fourth Planning and Development District 
is presently developing a regional solid waste management plan for its 
ten-county area including Spink County. The South Dakota Department 
of Environmental Protection is presently developing a state-wide 
solid waste management plan. This plan will identify the responr, ibi 1 itie 
and authorities for implementation of the State Plan. The plan will 
address the establishment of new open dumps, and the closing or upgrading 
of all existing open dumps, among other items. The U.S. Environment^' 
Protection Agency (USEPA) is required by law to inventory all disposal 
facilities or sites In the United States which are open dumps, and 
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those that do not meet the requirements of a sanitary landfill will bp 
placed on a compliance schedule for closure or upgrading. 

The problem is to determine where atd how many sites exist and 
ootential pollution problems. Many of the sites are near small towns 
but many more are in rural areas where generally a small group of people 
dump waste. The conventional way to find the sites would be a "wind- 
shield survey" of the area. Unfortunately this method is very costly, 
extremely time consuming and ineffective in locating sites not near 
main roads. Remote sensing techniques, which are generally cost 
effective for large area reconnaissance, were used to locate and map 
open and abandoned dump sites in Spink County. 

STUDY AREA 

Spink County encompasses approximately 968,000 acres. There are 
11 communities and, according to the 1970 census, a population of 
10,595. The triajor land use in this area is agricultural . The James 
River and the James River Basin, with numerous meandering streams, 
are partially ’ocated within Spink County. 


PROCEDURES 


Existing and newly acquired 1979 NASA high altitude color infrared 
photography was interpreted to ide; ’.ify dump sites in Spink County. 
Initial interpretation was pe>"formed using June 1975 1:60,000 scale 




67 


photography. The 1975 photography provided about 7555 coverage of the 
county. 

Additional interpretation was performed using the 1979 color infrared 
1:120,000 scale photography for the entire county. The photography was 
acquired in conjunction with other NASA projects. The quality of the 
data was considerably less than that of the 1975 data. Segments of 
many frames were overexposed and some clouds and cloud shadows were 
present. 

Ground truth information for known existing open dumps was provided 
by the Fourth Planning and Development District. The information 
provided was the location of the dump sites near each of the towns 
within the county. The ground truth information helped train an 
Interpreter to recognize possible dump sites. 

The location of possible dump sites was determined by using two 
key indicators. They were the presence of a mottled texture and tone, 
which was distinct from the usually homogeneous texture and tone of 
the surrounding area, and the presence of a road or trail used by 
vehicles to deliver the waste. Another indicator was the white tone 
of metal containers which was evident on the imagery. Both dates of 
imagery were interpreted using a 7x magnifying eyepiece. 

The accuracy of the interpretation was determined by field checking. 
Field checking was performed by trainee observers flying at low altitude 
in a small plane from site to site. The light plane method was chosen 
because of the difficulty in reaching sites not near roads and because 
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considerable time and money was saved In comparison to driving to each 
site. It took only three hours to verify 72 sites. 

Other resource Information concerning each site's pollution potential 
such as land use was also considered. The proximity to surface water 
was the most Important from a pollution standpoint since contamination 
of water supplies corld be very dangerous. In addition contamination of 
livestock forage by residual pesticide and herbicide chemicals may 
present a danger to the livestock. The basic land use surrounding 
each dump site was interpreted from the 1979 high altitude photography. 

The land use was grouped into four basic categories: cropland, pasture, 
riparian, and farmstead. The distance to surface water was calculated 
from US6S 7 1/2 minute topographic maps. The maps were used because 
the poor quality of the 1979 film made it difficult to accurately 
determine the drainage. The surface water bodies were either streams, 
lakes or marshes. No Information concerning location of shallow ground- 
water or any other geologic information was available. 

RESULTS 

A total of 146 possible dump sites were interpreted from the 1975 
and 1979 Imagery. A total of 72 sites were field checked, of which 33 
were actual garbage dump sites (Fig. 1). Garbage dumps were composed of 
cans, bottles, paper, and other household refuse. Another 20 sites were 
rubble piles composed of lumber, automobiles, and machinery (Fig. 2), 
and 19 were not dump sites of any kind. Only one dump site which had 
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not been interpreted from the imagery was found during the field 
checking. This indicates that the interpretation detected almost 
every probable dump site. 

Based on the interpreted dump sites, which were field checked, the 
level of accuracy for detecting probable garbage dump sites is 46, . 

If rubble sites which closely resemble garbage dump sites were 
included the accuracy increases to 75', 'i. 

A comparison of the accuracy results for dump sites interpreted 
from both scales of imagery was made to determine if one would be more 
advantageous. Of the 33 confirmed dump sites 19 were interpreted from 
the 1:60,000 scale imagery and 24 from the 1:120,000 scale photography. 
Consideration must be given to the fact that the 1:60,000 scale 
photography provided only 75% cove age of the county, the quality of 
the 1:120,000 scale imagery was poor and some sites could have been 
created after 1975. These three considerations make it difficult 
to make a complete evaluation of photographic scale. 

The dump site inventory was put in tabular form. The tabular 
information included a site number, the location based on a township 
and range system, the name of the USGS 7 1/2 minute quadrangle map 
in which the site is located, whether or not it was verified, the 
surrounding land use type, proximity to surface water, type of 
surface water and date of imagery from which the site was delineated 
(see Appendix A). The sites were also annotated on the appropriate 
USGS 7 1/2 minute topographic map to facilitate finding the sites while 
in the field. 
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SUMMARY AND CONCLUSIONS 

Interpretation of two dates and scales of imagery delineated 
146 potential open or abandoned dump sites. The inventory of sites 
and anrillary data were put into tabular form. Local qoverninent 
agencies were provided t‘ e Mata to a.ssist the refinement of the Spink 
County solid waste management plan and as input into the implementation 
of the State Plan. 

Despite the acceptance of the data as being accurate, timely and 
cost effective (see letter in Appendix B) the data have not been used 
extensively. Tne primary reason is the sensitive nature of solid waste 
management in Spink County. Spink County riicently withdrew from a 
multi-county solid waste program. Hence establishment of a new waste 
disposal site in Spink County has become a political issue. Hopefully, 
circumstances will change and the data will be used for planning purposes. 

Efforts will continue to expose those in the state and federal 
government to procedures developed. The study provides an economically 
and technically feasible method to inventory waste disposal sites. All 
indications are that similar inventories will be conducted nationally at 
some time in the future. 
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Mr. Jeff tidcnshink 
Remote Sensing Institute 
South Dakota St.Ue 
University 
Brookings, SO 57006 

Dear Jeff: 


The Fourlii Planning and Development District would like to express 
its appreciation to the Remote Sensing Institute at South Dakota State 
University Cor its efforts in the "Survey and Analysis of Potential 
Pollution fron Open and Abandoned Solid Waste Duiiip Sites Using Remote 
Sensing Techniques" in Spink County, South Dakota. 

The data compiled in the report appears to be most impressive in 
terms of accuracy, expedience, and cost effectiveness for a project of 
this nature. 

As to how this data will be utilized in Spink County is a matter 
that the local elected officials will have to examine and address. 

It would appear that the Reniote Sensing Institute, by the results 
of the Spink County Project, could take an active role related to 
Section 4005 of the Resource Conservation and Recovery Act (P.L. 94-580) 
in inventorying open and abandoned dump sites. 


If this office could be of any further assistance, please feel 
free to contact us. 
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1NFLUE^XE OF SOIL REFLECTANCE ON 
LANDSAT SIGNATURES OF CROPS 

INTRODUCTION 

The world economy in past decades to the present has exhibited 
an ever-increasing influence on modem American agriculture (Hueg, 1979). 
Agencies and individuals connected with agriculture now, more than ever, 
require a Voriety of information on a global basis. Accurate estimation 
of crop areas and yields are among the data which are crucial to decision 
makers at all levels. 

Estimation of crop area in the U.S. and numerous other countries 
has been accomplished largely by two major procedures: 1) the Area Frame 
concept developed by the Economics, Statistics, and Cooperatives Service 
(ESCS) which employs statistiCcHy-based ground data collection and 
expansion for agricultural data, e.g. livestock, crops, etc. (Von Steen 
and Wigton, 1976); and 2) interpretation and classification of Landsat 
data (Paarlberg et al., 1978; Malila and Gleason, 1977). The ESCS 
has been involved in both of the above procedures (Starbuck, 1977; 

Wigton and Bormann, 1978). 

The influence of soils, topography, and phenology has been 
cited as contributing to variability of crop signatures on Landsat 
data (Tucker and Miller, 1977; Wiegand et al., 1977: Westin and 
Lemme, 1978). Recent studies have investigated methods to negate the 
influence of the factors listed above by stratification (Dalsted 
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et al.i 1979; Padmanabhan . 1980) and by transformation of Landsat 
digital data (Kauth and Thomas, 1976). 

Holmes et a1. (1979) reported moisture stress and soil brightness 
as likely causes of altered spectral signatures of wheat on Landsat 
data. They also noted that stratification was a "useful tool" leading 
to "greater overall efficiency In producing estimates . . .". Schubert 
et a1. (1980) stratified five Landsat scenes In western Canada by 
Uniform Productivity Areas (UPAs), which were based on soil properties. 
In order to Improve sampling strategy. They reported a stratification 
technique that both aided In the selection of Landsat training data 
and related to the yields of spring wheat. The yield data within each 
UPA were based on 10-year production records. 

The objective of this study was to Investigate the statistical 
Implications of stratification upon crop signatures from Landsat CCT's. 
The stratification was based on manual photo-interpretation of multi - 
date Landsat coverage of an area In eastern South Dakota. 

STUDY AREA 

The region of Investigation Is shown In Fig. 1. Physiographic 
diversities and agricultural activities of this region have been 
detailed (Westin and Malo, 1978). Corn, wheat, oats, and barley are 
the major crops In this region. Sunflowers, however, are rapidly 
becoming a major crop In certain parts of this study area. 
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This region Is characterized by hot summers and co'd winters. 
Precipitation occurs predominantly In the spring and early simmer 
months. 


MATERIALS AND tCTHODS 

The study area was divided Into 10 strata by photo- Interpretation 
of 1:!>00,000 1978 Landsat prints (14 May, 3 August. 8 September 1978 
scenes). Locations of representative ground samples within each stratum 
were deten.iined by ESCS techniques (Von Steen and Wigton, 1976). The 
samples or segments were generally one-quarter section In area and 
totaled 255 for the entire region. The number of segments per stratum 
was weighted according to the areal extent of the respective stratum. 

Ground Investigation was accomplished by an on-site description 
of the land use at each of the segments; this data collection took place during 
the last week of July arid first week of August, 1979. Aerial photography 
(1:15,840 scale) was used as a base for outlining field boundaries and 
notating crops and crop phenology (Tabic 1) at each segment. The aerial 
photography used In the field Investigation varied In date, the majority 
being In the 1975 to 1979 period. 

The August 25 data were chosen for the Initial analysis because 
of Interest shown by ESCS personnel In using these data along with the 
RSI-acquired field Information for verification of their sunflower 
acreage estimates. Time limitations precluded the use of the July 20 
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Table 1. Crop phenology (ad hcc) Is listed by code and characterizing 
features. These stages were noted at the time of ground data 
collection, approximately two to three weeks before the 
25 August Landsat data. 


Code 

6 

7 

8 
9 

10 

11 

12 


Stage (row crops) 

.3 - .6 m In height 
.6 • 1 .2 m In height 
1 .2 • 1 .8 m in height 
> 1 .8 m In height 
0 - 25X bloom or tassel ling 
25 • 50S bloom or tassel ling 
> 50* hloom or tassel ling 


data. The 1ate*season date of the Landsat data utilized limited this 
study to the analysis of sunflower and com data. 

The segment data were extracted from the Landsat CCT using the 
ESCS-deve loped EDITOR computer software package (Steinberg, 1979*, 
Institute for Advanced Computation, 1978), the Bolt, Baranck and 
Newman, Inc. DEC PDP>10 computer, and established ESCS techniques for 
map digitization and data verification. A master data set was created 
which contained the pixel values for MSS bands 4 through 7 and tagged 
each pixel by crop, strata, segment, field and stage of growth. Sub> 
sets of this master data set were created for the com and sunflower 
data. 

Various statistical procedures from the Statistical Analysis 
System (SAS) software package (SAS Institute, 1979) were used to analyze 
the com and sunflower data sets. These procedures Included MEANS 
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(means and other descriptive statistics). GLM (general linear tiiodels - 
analysis of variance), and DUNCAN (Duncan's multiple raiiqe test). A 
nested model was used for the ULM analysis with tiie segnieni and field 
variables being random and the strata and stage variables being fixed. 

RESULTS AND DISCUSSION 

The approach of this study was designed to develop a Landsat 
stratification procedure for application to other crop>producing regions 
of the world. To reduce bias (i.e. prior knowledge of the soils in the 
region) only Landsat imagery was used to establish the strata. An 
operational approach v/ould utilize Landsat and physiography data in 
concert. 

The strata boundaries were established using three Landsat dates; 
the August and September dates gave an indication of land use and land 
capability while the May date showed soil tones (Fig. 2). The May date 
was given priority for boundary location followed in order by the 
September date and August date. 

Since an August 25. 1979 CCT was used^ jnly corn and sunflowers 
were considered in the analysis. Analysis and interpretations in 
conjunction with CSCS are still in progress. 

Com 

Mean digital count data by strata are shown in Table 2. All but 
stratum 1 are well represented in number of observations. Standard 
deviations values within each band were consistent with the exception 
of strata 6 and 9 for MSS 6 and 7; this indicated a reasonable level of 
accuracy from the field identifications and data extractions. 





Table 2 . Landsat digital means (25 August 1979 CCT) are listed by strata 
and Landsat band for corn. 


MSS 4 

Strata 


Mean 

Standard Deviation 



1T7 

Tff.T 

D'.T 


2 

312 

18.5 

2.4 


3 

641 

17.1 

1.3 


4 

617 

16.9 

1.3 


5 

532 

16.1 

1.3 


6 

342 

17.2 

1.9 


7 

1321 

15.0 

1.6 


8 

449 

14.5 

1.2 


9 

357 

15.1 

1.5 


10 

366 

15.2 

1.2 

MSS 5 

1 

117 

17.7 

1.6 


2 

312 

18.2 

4.5 


3 

641 

16.4 

2.4 


4 

617 

16.4 

2.1 


5 

532 

15.5 

2.0 


6 

342 

16.4 

2.4 


7 

1331 

14.3 

2.8 


8 

449 

13.6 

1.5 


9 

357 

14.9 

3.1 


10 

366 

14.3 

1.8 

MSS 6 

1 

117 

43.0 

3.9 


2 

312 

46.9 

4.2 


3 

641 

43.0 

3.1 


4 

617 

44.9 

2.9 


5 

532 

43.7 

5.6 


6 

342 

48.4 

12.6 


7 

1331 

42.5 

4.8 


8 

449 

39.6 

4.8 


9 

357 

39.9 

7.8 


10 

366 

42.1 

4.1 

MSS 7 

1 

117 

43.8 

5.4 


2 

312 

46.9 

4.2 


3 

641 

44.8 

4.0 


4 

617 

47.7 

3.8 


5 

532 

45.4 

7.1 


6 

342 

52.4 

16.0 


7 

1331 

45.9 

6.1 


8 

449 

42.6 

6.1 


9 

357 

42.1 

10.3 


10 

366 

43.9 

5.6 


^Observations, Landsat pixels 
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The mean digital data were also subdivided by crop stage (Table 3). 
This subdivision, which had 5 crop stages, did not recognize strata and 
provided a larger sample size/unit than did the subdivision by strata. 
Several general trends are visible in the data. The more mature classes 
(i.e. units 9 and 10) had higher digital values in Landsat band 4 and 
band !>. In Landsat band 6 and band 7 unit 10 had the lowest mean digital 
value, and unit 9 had the highest mean value. While the ground data 
recorded crop stage approximately two to three weeks in advance of the 
25 August CCT, the data followed predictable trends. The greater the 
crop maturity the higher the digital count in the visible spectrum and 
the most advanced maturity (i.e. unit 10) had the lowest digital count 
in the near-infrared spectrum. 

Table 3. H^an digital counts of corn data are shown by crop stage 
and Landsat band. 


Stage 

N 

MSS 4 

MSS 5 

MSS 6 

MSS 7 

6 

340 

15.2 

13.9 

43.8 

47.0 

7 

434 

15.7 

14.7 

43.6 

46.7 

8 

1069 

15.6 

14.9 

42.9 

45.7 

9 

518 

16.1 

15.4 

44.9 

47.2 

10 

2624 

16.3 

15.9 

42.6 

45.1 


Analysis of Variance (ANOVA) results Indicated significant 
differences among variables (Table 4). Since one or more stage and 
stratum were indicated as significantly different among the variables, 
a multiple comparison test was utilized to locate those differences by 
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Table 4. AfJOVA results are shown for corn from the August 25, 1979 
CCT. A fixed model was employed. 




Independent Variables 


Dependent 

Variable 

Stage(5)''’ 

F Values 
Strata(9) 

Stage* strata(19) 


Band 4 

20.68** 

149.20** 

33.55** 


Band 5 

28.31** 

45.62** 

32.26** 


Band 6 

81.61** 

163.36** 

147.75** 


Band 7 

68.38** 

145.97** 

143.47** 



t 


Degrees of freedom 


strata (Table 5). Strata are divided Into five to seven significant 
groups, dependent upon Landsat band. The Interplay of the crop, soil, 
crop stage, and wavelength sensitivity combined the various strata In 
different configurations for each MSS band. Further statistical testing 
Is needed to establish one system of strata separation/combination, e.g. 
multivariate analysis, covariance analysis, clustering, etc. 

Duncan's New Multiple Range Test (DMR) was also used to Indicate 
the significant groupings of digital count means for corn stage (Table 6). 
Landsat bands were shown to have differing sensitivities to crop stage. 

The advanced maturity of stage 10 resulted in high digital count in 
the visible bands (i.e. 4 and 5) and a low count in the near-infrared 
band (i.e. 6 and 7). Since each crop stage did not occur within each 
stratum, the data set Is Incomplete. Further multiple comparisons are 
needed to examine the Interaction of strata and crop stages. Phenology 
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Table 5. Duncan's New Multiple Range test for August 25 corn data by 
strata. Mean digital count decreases in magnitude from top 
to bottom. 


MSS 4 

Att O 


MSS 5 


MSS 6 


MSS 7 


B 1 


Observations are aligned with the order of strata in band 4, i.e 
unit 2 - 312 pixels, etc. 

Units with the same letter are not significantly different (.05 
probability level). 
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Table 6. Duncan's New Multiple Range test for August 25 corn data by 

crop stage. Occurrence of stage within strata Is also listed. 


Band 4 

Band 

_5 

Band 6 

Band 

_7 

A* 

A 

10 

A 

9 

A 

9 


A 




A 


B 9 

A 

9 

B 

6 

B A 

6 




B 




C 7 

B 

8 

B 

7 

B 

7 

C 

8 


B 




C 8 

B 

7 

B 

8 

C 

8 




B 




D 6 

C 

6 

B 

10 

0 

10 


^Unlts with the same letter are not significantly different (.05 level 
of probability). 

^^Occurrence of crop stage by strata Is as follows: stage 6 - 5,7,9,10; 
stage 7 - 2,4,6,9,10; stage 8 - 3,4,5,6,7,8,9,10; stage 9 - 1,2, 3, 7, 8 
9,10; and stage 10 - 1 ,2, 3, 4, 5, 6, 7, 8. 


differences within one -p across a region cause classification problems. 
Further strata comparisons with phenology held constant may provide 
additional information about the strata/crop stace interrelationship. 

Sunflowers 

Mean digital counts of sunflowers by strata are shown In Table 7. 
Strata 1, 5, and 10 were not represented. Several strata had a low 
number of observations; consequently, sensitivity of statistical tests 
will be reduced. However, the standard deviations appeared to be 
reasonably low. 

Mean digital counts of sunflowers by crop stage are shown in Table 8. 
The advanced maturity of stage 12 was apparent by the low digital values. 
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Table 7. Landsat digital means (25 August 1979 CCT) are listed by strata 


ta ndsat band for sunf lowers . 


K5 S 4 

Strata 


Mean 

Standard Deviation 


2 

763 

20.0 

1.3 


3 

10 

18.3 

0.7 


4 

73 

18.1 

1.2 


6 

113 

19.0 

1.1 


7 

357 

17.3 

1.2 


8 

51 

17.1 

0.7 


9 

95 

17.6 

1.3 

HSS 5 

2 

763 

19.1 

1.8 


3 

10 

17.1 

1.0 


4 

73 

18.3 

2.1 


6 

113 

18.2 

1.2 


7 

357 

16.3 

1.8 


8 

51 

16.2 

1.1 


9 

95 

17.4 

1.4 

MSS 6 

2 

763 

63.1 

5.9 


3 

10 

59.6 

4.2 


4 

73 

57.6 

5.5 


6 

113 

57.5 

6.9 


7 

357 

53.4 

9.0 


8 

51 

60.1 

6.0 


9 

95 

52.7 

8.4 

MSS 7 

2 

763 

63.5 

6.5 


3 

10 

63.0 

5.2 


4 

73 

60.1 

6.2 


6 

113 

58.1 

7.4 


7 

357 

54.7 

10.1 


8 

51 

62.4 

5.6 


9 

95 

51.8 

8.8 


^Observations, Landsat pixels 
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Table 8. Landsat digital means (25 August 1979 CCT) are listed by 
crop stage and Landsat band for sunflowers. 


Stage 

N 

f4SS 4 

MSS 5 

MSS 6 

14SS 7 

7 

196 

18.6 

19.0 

52.6 

52.2 

8 

344 

19.5 

18.4 

62.4 

63.1 

9 

450 

19.4 

18.3 

62.8 

63.9 

10 

- 

- 

- 

- 

- 

11 

330 

18. C 

17.6 

59.0 

59.9 

12 

142 

17.2 

16.0 

50.0 

51.0 


Analysis of variance and multiple comparison tests are In progress. 
A comparison of these results with the results from the corn data will 
provide an additional basis for evaluation of the procedure. 

The sunflower ground data proved to be very valuable to the FSCS. 
They were able to expand the Information to produce an area estimate of 
sunflowers, a crop whose rapid growth In the area was of significant 
Importance. 


SUMflARY AND CONCLUSIONS 

A stratification procedure for partitioning Landsat data Into 
homogeneous regions relative to crop identification was Investigated. 
Strata were evaluated on the basis of corn and sunflower data from an 
August 25, 1979 CCT; the study area was a six-county area In eastern 
South Dakota. 
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Results are not yet complete, but the following conclusions can be 
made relative to com: 1) utilizing a fixed ANOVA model, there are 
differences among strata and crop stage In the digital count means; 

2) five to seven strata, dependent upon Landsat band, were recognized 
as significantly different; 3) two to four groupings of crop stages out 
of a possible total of five crop stages were determined by multiple 
comparison. This Indicated differing sensitivity of Landsat bands to 
crop stage. Further analysis Is underway, I.e. holding crop stage 
constant while testing strata. 

Analysis of sunflower data Is In progress. The results will be 
compared to the corn results. Sunflower ground Information was expanded 
to provide the ESCS with an area estimate of this crop which has rapidly 
Increased In Importance In South Dakota. 
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APPLICATION OF REMOTE SENSING IN 
THE NATIONAL MODEL IMPLEMENTATION PROGRAM, 

LAKE HERMAN WATERSHED 

INTRODUCTION 

In 1978 Lake Herman watershed near Madison. South Dakota was 
selected as one of seven water resources systems In the United States 
to participate In the national Model Implementation Program (MIP). MIP 
Is a three year pilot cooperative program by the United States Department 
of Agriculture (USDA) and the United States Environmental Protection 
Agency (USEPA). The program Is designed to test and evaluate various 
water resources quality Improvement methods. 

Several South Dakota groups are cooperating In the Lake Herman 
MIP program. They Include: Lake Hennan Develop-ient Association, First 
Planning and Development District, East Dakota Conservancy District, 

South Dakota Department of Natural Resource Development, South Dakota 
State University Extension Service, USDA, Agricultural Stabilization 
and Conservation Service, USDA Soil Conservation Service (SCS), and 
the Remote Sensing Institute (RSI). 

At the outset of the MIP, remote sensing techniques were employed 
as a means of monitoring land use and hydrologic changes throughout 
the watershed. In 1978, NASA and SCS funded data collection and 
analysis and provided technical expertise for establishment of an 
Inventory of the baseline condition of the watershed. SCS Is the 
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major user of the data. Baseline landcover viata for the 1978 <jrow1n9 
season and detailed soil survey data were coded Into RSI's Area 
^source ^alysls System (AREAS), a computerized cellular Information 
system (Wehde, 1979). The data were analyzed to produce tabular and 
map Information concerning location, quantity and severity of soil 
erosion (Myers, et al. 1979). 

The results of the first year of activity began to address several 
of the MIP objectives. Among them are the need to establish the 
baseline condition of the watershed, Identify land In need of treatment, 
Identify sites for sediment control structures, monitor changes within 
the watershed resulting from MIP and to document the achievements of 
MIP. Continual updating of the landcovcr data Is necessary to address 
many of the objectives of MIP. The results of the second year of 
activity continued to address these objectives. 

PROCEDURES 

Initially, It was intended to utilize NASA high altitude photography 
acquired during the 1979 growing season for Interpreting the landcover. 
Unfortunately, weather conditions prohibited data collection during 
the time period requested. Consequently, 70 nm 1:61,000 scale color 
and color Infrared photography was acquired on 29 August using the RSI 
low altitude aircraft. 

Groundtruth data for the watershed were collected prior to the 
Interpretation of the photography. Twelve landcover categories 
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Identical to those used for the 1978 data were Interpreted. The 
groundtruth data were obtained by recording the landcover type that 
occurred in each quadrant at most road intersections within the water- 
shed. The Information was used to train the interpreters and verify 
interpretation results. 

The 1979 landcover data were computerized using a method similar 
to the 1978 process (Myers et a1. 1979). The photographic data were 
enlarged to approximately 1:11,000 using an International Imaging 
System (I S) color additive viewer.* The enlarged data were overlaid 
with a cell grid. Each cell represented .25 hectares (.625 acres). 
Each cell was classified, based on category dominance within the cell. 
The digitized data were verified for all category boundaries and 
classes. 

Lan d Treat ment 

Several hundred acres of land were terraced during the first 
year of MIP. Much of the land was identified as in need of treatment 
from maps generated from analysis of the 1978 landcover data and 
detailed soil data within the Information system. The land treatment 
data set was updated to reflect the improvements. Infuntiation 
concerning the location of the newly constructed terraces was provided 
by SCS. A revised map of land in need of treatment was produced for 
SCS. 


•Inclusion in this report of registered tradenames or trademarks does 
not constitute an endorsement by the author or the Remote Sensing 
Institute. 
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Sediment Del Ivery Modeling 

The latest appltcation of the infonnat Ion system is l>c*imj <lircctc<J 
toward tiiodeling sodinwnt delivery to the lake. SCS has dellrtcateU four 
major tributary areas which are delivering sediment to the lake. The 
subdivision of the natural surface drainage network of streams accounts 
for three areas. The fourth area is the innwdiate lake shore. The 
fourth area and associated drainage network was delineated and encoded 
Into the information system. 

The tributary areas, drainage network, land cover and capability 
class data sets were used in the two approaches which wore investigated 
to estimate sediment delivery to the lake. The first approach was to 
assign a weight to each cell of cropland on capability class HE, 1 1 IE, 
and IVE within the watershed. Cropland on these classes is considered 
to be the primary source of sediment by SCS. Tiie weight was based on 
the shortest distance of the cell to the lake shore. Consequently, a 
cell on the lake shore would be weighed much more iieavily than one 
on the outer extremity of the watershed. Theoretically the potential 
for sedimeit reaching the lake would decrease as the distance to the 
lake increased ( Fig. 1). 

The estimated sediment yield rate for the soils within tne 
watershed was estimated in tons per acre. The estimates were made 
by SCS using the rainfall factor (R), erodibility factor (K), slope 
length and gradieni factors (LS), vegetative cover factor (C), and 
the erosion control practice factor (P) of the universal soil loss 
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equation, estimates were calculated for row crops, small grains, and 
grassland on capability class HE, HIE, and IVE land. The estimates 
do not take Into consideration any existing land treatment. 

The SCS estimates of sediment yield rate were used In a formula 
to calculate total sediment delivery. The formula Is: 

Wlu \c/lu - T/A 

where: 

\c/lu Is the acreage of a landcover class on a capability class 

YR Is the yield rate 

Wp Is the weight based on proximity to the lake 

The second approach Involved weighing each cell of cropland and 
grassland. This time however the cells were weighed on the distance 
to drainage rather than the lake Itself. The basic assumption In the 
model was that any sediment entering the drainage would eventually 
enter the lake (fig. 2 ). 


RESULTS 

Sediment Control Structure Site Planning 

One of the major MIP activities of SCS is the planning and 
construction of sediment control structures. At least three major 
structures and several minor structures will be constructed. During 
the summer of ’979 SCS was involved in planning a major sediment 
control structur»» in the drainage entering the southern tip of Lake 


Herman. 
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The SCS utilized specific map and tabular information derived from 
AREAS throughout the planning and design of the sediment control structure. 
One of the major uses of the information was to determine the necessary 
size of the structure. The SCS utilized landcover, land tret, ted, land 
in need of treatment and soil capability class information to estimate 
potenti'*. runoff and soil erosion. 

Another use of the information was to determine potential sites. 

Sites were determined from a slope map of the watershed generated from 
the detailed soil data. The map delineated areas of steep slopes suitable 
for construction sites. The final sites were selected from field surveys. 
Figure 3 delineates the areas of steep slopes with existing and proposed 
sites annotated. 

Sediment Delivery Modeling 

Two sedimentation models have been evaluated to date. The first 
approach which was based jn the distance of a parcel of land from the 
lake shore, utilized the SCS sediment yield estimates for cropland 
and capability classes (Table 1) to estimate sediment delivery to the 
lake. 

Table 1. Estimated sediment yield (tons/acre) for landcover type and 
capability class. 



HE 

IIIE 

IVE 

Row Crops 

5.3 

10.2 

15.3 

Small Grains 

1.8 

3.9 

6.3 


Grassland 


.2 


.9 


1.8 
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The estimate of sediment delivery to the lake from the first approach 
for all four tributary areas was 42,800 tons or 32.7 acre feet. This 
estimate was relatively close to the results obtained from a study 
conducted by USDA Science and Education Administration. Sedimentation 
Laboratory (ESCS, 1979). The study involvea collecting several core 
samples of sediment from the lake bottom. The age and rate of 
sedimentation was determined by measuring the amount of Cesium 137 in 
the core samples. The presence of the radioactive Cesium 137 was the 
result of nuclear testing conducted in 1954. The study measured a 
25 year average sediment delivery rate of 22 acre feet annually. The 
estimate of the first approach would be closer to the ESCS measurement 
if existing land treatment is considered. 

A breakdown by tributary area indicated as expected that the 
largest area was the greatest contributor. However, the breakdown 
showed that the average yield in tons per acre from 'ropland and 
grassland was the greatest in the second largest area. This indicated 
to SCS th'.t perhaps a more intensive soil conservation effort was 
needed in that area. 

The second approach which was based on the distance to drainage 
approach derived a sediment yield estimate of 72,900 on 55.8 acre feet. 
This is nearly 58 percent higher than the first estinate. The results 
of the second approach lead one to believe that the assumption that 
all sediment entering drainage reaches the lake is incorrect. Con- 
sequently, a new approach is being investigated which not only weighs 


the distance to drainage but also weighs the distance of the drainage 
point from the lake along the drainage path. It Is anticipated that 
this method will provide a more reasonable estimate. 

CONCLUSIONS 

Over the past two years SCS has been using a geographical 
Information system to address the objectives of MIP. The objectives 
Include establishing the baseline condition of the watershed. 
Identifying land In need of treatment, constructing sediment control 
structures and documenting the achievements of MIP. 

The map and tabular Infonnatlon derived from the Information 
system have greatly facilitated the aspects of MIP. SCS has been 
able to utilize the map and tabular data to reduce the amount of field 
work and concentrate In areas where land treatment Is most necessary 
and to delineate areas most suitable for sediment control structures. 

The sediment yield model has provided new insight Into the effects 
of land management on sediment delivery. Although It Is difficult 
to assess the accuracy of the model repeated use of the same model over 
time will provide a relative evaluation of the overall reduction of 
sediment yield as a result of land treatment or other management 
practices. Further refinement will Improve the usefulness of the 
model for overall watershed management and sediment control structure 
es 
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SCS has realized the benefits of utilizing natural resource 
data acquired from remote sensing and other techniques In an 
Information system approach. Consequently, SCS Is providing major 
funding support to continue the work through the remaining year of MIP. 

In the future the landcover and land treatment data sets will be 
updated to facilitate the evaluation of the Impact of the MIP program. 
Crop rotation and land treatment practices will be studied and further 
Improvements will be made within the watershed. In general. It Is 
anticipated that the use of natural resource Information systems of 
this design could be used to address land and water quality management 
problems nationwide. 
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SIX-MILE CREEK IMVESTIGAT lOH FOUOW-ON 

f • 

INTRODUCTION 

An invcsti(j<it1un was initiated in I97H t'l evaluate the geo- 
hydroloqy and the on/iromnent of Six-Mile Creek Watershed near 
Brookings, South Dakota to evaluate the impact of a proposed dam and 
reservoir site (Myers, et al., 1979). Thermal and color IR data were 
used by SCS and RSI personnel to evaluate the ground water situation, 
complete a land use inventory, and characterize erosion potential for 
preparation of an environmental impact statement. Specific objects 
were to: 

1) Select sites for drilling observation wells 

2) Evaluate water-table depths 

3) Evaluate soil noisture 

4) Evaluate erosion potential 

These objectives were met in FY 79. An extension of the Six-flile 
Creek Investigation was requested so that results of the drilling 
program (objective 1) could be fully documented. Refer to the 1979 
Annual Report SDSU-RS 1-79-1 4 for a complete discussion of the 
Investigation. 


RESULTS 

A ground water investigation of the proposed dam and recreation 


site o.j Six-Mile Creek Watershed (Fig. 1) is required to evaluate 



potential seepage losses. The Investigation consists of drilling and 
logging of observation wells to evaluate the vertical extent and layering 
of permeable aquifer material. Wells are Installed so that ground water 
heads can he measured Independently In each aquifer at multiple lateral 
locations. Relative and actual permeability and boundaries of 
permeable material must be known to design seepage cutoff and drainage 
control under the structure. 

The geohydrologic Investigation of the site was Initiated In 197U 
with the drilling of 11 observation wells along the proposed structure 
site and near the northern end of the pool area. Potential sites for 
drilling remaining wells were selected on the basis of the remote 
sensing imagery acquired in early September, 1978 (Myers et al., 1979). 
Thermal inagery of the recreation site revealed a broad cool pattern 
surrounding Six-Mile Creek (Fig. 2). The shape and location of the 
apparent thermal anomaly were typical of thermal patterns produced by 
shallow aquifers. 

Test holes 805, 806, and 807 (Fig. 2), later renumbered 801, 802, 
and 803, respectively by SCS, were drilled in August 1979. Results 
of the drilling are summarized in Table 1 and Fig. 3-5. Drilling 
verified the presence of permeable alluvium at all three locations. 

Depth to water ranged from 1 to 5 feet below the soil surface. 

Observation wells and piezometers were installed at the three 
locations to monitor ground water fluctuations and head pressure. 
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SUMflARY 


This completes the S1x«M11c Creek Investl9at1on. The attached 
letter from the State Conservationist suninarlzes the cooperative study 
between the Soil Conservation Service and the Remote Sensing Institute 
This project hes been a successful demons traticn of the use of remote 
sensing In watershed planning and development. Information obtained 
from the Investigation will be used to prepare the Preliminary 
Investigation Report (PIR) which will form the basis for the final 
Working Plan and Environmental Impact Statement. The PIR will be 
submitted to NASA upon completion. 
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Mr. Victor layers • 01 rector 
Roniott Sensing Institute 
South Dakota State University 
brookings* South Dakota S/007 

Dear Mr. Myers: 

Me wish to thank you and your staff for the fine cooperotion and help 
you have given us. 

The cooperative study between the flriiute Sensing Institute and the Soil 
Conservation Service on Six Mile Creek WatersheJ has been extremely 
useful and timely in the preparation of the watershed plan for flood- 
damage reduction and the development of water-based recreation. 

Identification of land use and land-treatment needs arc integral parts 
of the environmental assessment of the watershed. The remotely sensed 
data obtained by RSI is being used in this idenli ficatiun process. 
Digitized Imagery has been analyzed and put in tabular fomi by RSI 
for soils mapping units and land used. The land use data has been 
superposed on watershed maps by RSI. SCS will pick out key soils which 
RSI will superpose on the land use maps. This data will be used by SCS 
to locate and quantify flood-prone areas and land treatment needs. 

The proposed multipurpose reservoir designated SM-IB, located north of 
White, South Dakota, is being investigated as both a flood-water 
retarding and recreation reservoir. Those combined purposes require 
both safety and watei holding capability considerations. Understanding 
the ground-water regime of the watershed and present and potential 
subsurface-water movement to and from the proposed reservoir are 
required in both safety and water holding capability considerations. 

Infrared and near surface temperature scan data obtained by RSI were 
used to locate the boundaries of near surface permeable alluvium. 
Preliminary drilling was used to identify the subsurface materials. 

The location of the drill holes was selected using both the remotely 
sensed data and traditional site exploration methods. Ground-water 
observation wells and piezometers v>ere installed to monitor the 
fluctuations in ground-water levels and head pressure. The monitoring 
is being continued by SCS to obtain adequate data on seasonal changes. 
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The renotely sensed dJt<i has made It possible to reduce the time end 
expense of subsurface drilliioj needed to identify the earth material 
characteristics. SCS has been able to more efficiently study the ground 
water conditions because the remotely sensed data has boin used in 
decisions on the location aid nunter of ground-water observation 
stations. 

Cost and t1n»e efficiency luive been ini(>f‘Oved by the use of remotely sensed 
data In the environmental assessment through reduction of time required 
for land use and land treatment inventories. 



R. 0. Swenson 
State Conservationist 
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cc: Kenneth Huber. WRPSL, SCS, Huron 

James Hyland. Geologisli SCS, Huron 
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